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Abstract 

This review aims at a critical discussion of the interplay between effective interactions derived from various many- 
body approaches and spectroscopic data extracted from large scale shell-model studies. To achieve this, our many-body 
scheme starts with the free nucleon-nucleon (NN) interaction, typically modelled on various meson exchanges. The NN 
interaction is in turn renormalized in order to derive an effective medium dependent interaction. The latter is in turn 
used in shell-model calculations of selected nuclei. We also describe how to sum up the parquet class of diagrams and 
present initial uses of the effective interactions in coupled cluster many-body theory. 
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1 Introduction 



Traditional shell-model studies have recently received a renewed interest through large scale shell model calculations 
in both no-core calculations ^1211310 for light nuclei, the IsOd shell the IpO/ shell and the 2sldQg-j/2 shell with 
the inclusion of the O/111/2 intruder state as well. It is now therefore fully possible to perform large-scale shell-model 
investigations and study the excitation spectra for systems with many basis states. With recent advances in Monte Carlo 
methods, see for example Refs. [S1[7II51I^ I1(J|[TT| |12 | . one is also able to enlarge the dimensionality of the systems under 
study considerably, and important information on e.g., ground state properties has thereby been obtained. 

An important feature of such large scale calculations is that it allows one to probe the underlying many-body physics 
in a hitherto unprecedented way. The crucial starting point in all such shell-model calculations is the derivation of an 
effective interaction, be it either an approach based on a microscopic theory starting from the free nucleon-nucleon (NN) 
interaction or a more phenomenologically determined interaction. In shell- model studies of e.g., the Sn isotopes, one 
may have up to 31 valence particles or holes interacting via e.g., an effective two-body interaction. The results of such 
calculations can therefore yield, when compared with the available body of experimental data, critical inputs to the 
underlying theory of the effective interaction. 

Clearly, although the NN interaction is of short but finite range, with typical interparticle distances of the order of 
1 '--^ 2 fm, there are indications from both studies of few-body systems and infinite nuclear matter, that at least three- 
body interactions, both real and effective ones, may be of importance. Thus, with many valence nucleons present, such 
large-scale shell- model calculations may tell us how well e.g., an effective interaction which only includes two-body terms 
does in reproducing properties such as excitation spectra and binding energies. 

This work deals therefore with various ways of deriving the effective interaction or effective operator needed in shell- 
model calculations, starting from the free NN interaction. Normally, the problem of deriving such effective operators and 
interactions are solved in a limited space, the so-called model space, which is a subspace of the full Hilbert space. The 
effective operator and interaction theory is then introduced in order to systematically take into account contributions 
from the complement (the excluded space) of the chosen model space. Several formulations for such expansions of 
effective operators and interactions exit in the literature, following time-dependent or time-independent perturbation 
theory 1141 115L I16L 1171 [TS| . The so-called Q-box method with the folded-diagram formulation of Kuo and co-workers 
[151 116j has been extensively applied to systems in nuclear physics, especially for the derivation of effective interactions 
for the nuclear shell-model. This method is however limited by the fact that the Q-box includes diagrams to a certain 
order in perturbation only. To go beyond third order in perturbation theory is very hard and there is no indications that 
the expansion converges. Thus, formulations like the coupled-cluster method or exponential ansatz llTllTfl. 20, 21 . 22, .2^ 
1251 12()l ITtI 1281 1^ 1301 1311 1321 \',VM IH^ 1351 llffij . which is much favored in quantum chemistry, the summation of the parquet 
class of diagrams ISZIOHIEnilinilllllHESlEl, although with few applications, and FHNC theory [111113113, offer the 
possibility of summing much larger classes of many-body terms. 

Ab initio methods like variational and diffusional Monte Carlo approaches [HI |H1 IIOL 1111 [T^ and no-core-shell model 
calculations ^ EJ El ^ for light nuclei have also been much favored in nuclear theory. 

In this work we will focus on the above-mentioned Q-box approach combined with the nuclear shell-model, the 
summation of the so-called parquet diagrams and Coupled cluster theory as three possible ways of studying the nuclear 
many-body problem. The Q-box has been introduced in Rayleigh-Schrodinger perturbation theory as the definition of all 
non-folded diagrams to a given order in the expansion parameter, in nuclear physics the so-called G-matrix. The G-matrix 
renders the free NN interaction V, which is repulsive at small internucleon distances, tractable for a perturbative analysis 
through the summation of ladders diagrams to infinite order. Stated differently, the G-matrix, through the solution of 
the Bethe-Brueckner-Goldstone equation, accounts for the short-range correlations involving high-lying states. Folded 
diagrams are a class of diagrams which arise due to the removal of the dependence of the exact model-space energy in 
the Brillouin-Wigner perturbation expansion. Through the Q-box formulation and its derivatives, this set of diagrams 
can easily be summed up. 

In addition to the evaluation of folded diagrams and the inclusion of ladder diagrams to infinite order included in 
the G-matrix, there are other classes of diagrams which can be summed up. These take into account the effect of long- 
range correlations involving low-energy excitations. A frequently applied formalism is the Tamm-Dancoff (TDA) or the 
random-phase (RPA) approximations. In their traditional formulation one allows for the summation of all particle- hole 
excitations, both forward-going and backward going insertions. This set of diagrams, as formulated by Kirson |47| and 
reviewed in Ref. .4^, should account for correlations arising from collective particle-hole correlations. Another possibility 
is to include any number of particle-particle and hole-hole correlations in the screening of particle-hole correlations. The 
inclusion of these correlations is conventionally labelled particle-particle (pp) RPA. It has been used both in nuclear 
matter studies 49, 5X1., ,51., .52^ and in evaluations of ground state properties of closed-shell nuclei |5;3i i54i i55j. Ellis, 
Mavromatis and Miither |54[ have extended the pp RPA to include the particle-hole (ph) RPA, though screening of 
two-particle-one-hole (2plh) and two-hole-one-particle (2hlp) vertices was not included. The latter works can be viewed 
as a step towards the full summation of the parquet class of diagrams. The summation of the parquet diagrams entails a 
self-consistent summation of both particle-particle and hole-hole ladder diagrams and particle-hole diagrams. Practical 
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solutions to this many-body scheme for finite nuclei will be discussed here. 

The coupled-cluster method originated in nuclear physics over forty years ago when Coester and Kiimmel proposed an 
exponential ansatz to describe correlations within a nucleus I2()j . This ansatz has been well justified for many-body 
problems using a formalism in which the cluster functions are constructed by cluster operators acting on a reference 
determinant 24 . Early applications to finite nuclei were described in Ref. |3T]. From that time to this, a systematic 
development and implementation of this interesting many-body theory in nuclear physics applications has been only 
sporadic. The view from computational quantum chemistry is quite different. In fact, coupled-cluster methods applied to 
computational chemistry enjoy tremendous success |21l 1221 051 126[I27| over a broad class of chemistry problems related to 
chemical and molecular structure and chemical reactions. The method is fully microscopic and is capable of systematic 
and hierarchical improvements. Indeed, when one expands the cluster operator in coupled-cluster theory to all A particles 
in the system, one exactly produces the fully-correlated many-body wave function of the system. The only input that the 
method requires is the nucleon-nucleon interaction. The method may also be extended to higher-order interactions such as 
the three-nucleon interaction. Second, the method is size extensive which means that only linked diagrams appear in the 
computation of the energy (the expectation value of the Hamiltonian) and amplitude equations. Third, coupled-cluster 
theory is also size consistent which means that the energy of two non-interacting fragments computed separately is the 
same as that computed for both fragments simultaneously. In chemistry, where the study of reactions is quite important, 
this is a crucial property not available in the interacting shell model (named configuration interaction in chemistry). 

A fourth interesting point involves a comparison of coupled-cluster theory and many-body theory. The computed 
energy using the coupled-cluster formalism includes a very large class of many-body perturbation theory diagrams. In 
standard many-body perturbation theory, one typically sums all diagrams order by order. The coupled-cluster approach 
essentially iterates diagrams so that one may discuss it in terms of an infinite summation of particular classes of diagrams. 
The infinite resummation means that the coupled-cluster theory is nonperturbative. In fact, the coupled-cluster energy at 
the single and double excitation level contains contributions identical to those of second order and third order many-body 
perturbation theory, but lacks triple excitation contributions necessary to complete fourth-order many-body perturbation 
theory; see e.g., the review article of Bartlett |21| . Later in this review, we will compare second and third-order many-body 
perturbation theory to coupled-cluster calculations. 

This work falls in eight sections. In the next section we present various definitions pertinent to the determination of 
effective interactions, with an emphasis on perturbative methods. The resummation of the ladder type of diagrams is 
then presented in sectional In that section we also discuss the summation of so-called folded diagrams which arise in the 
evaluation of valence space effective interactions. Further perturbative corrections are also discussed. Selected results for 
light nuclei in the IsOd and IpOf shells and for several medium heavy nuclei in the mass region A — 100 to A = 132 are 
presented in the subsequent section. 

The summation of the TDA and RPA class of diagrams and the so-called parquet diagrams is discussed in section 
whereas sectional is devoted to a discussion of the coupled cluster method. Section [3 presents a critical discussion of 
three-body effects in nuclear structure. 

We conclude in section |H1 with a discussion on the extension of the methods discussed in sections and El to weakly 
bound nuclei. 



2 Many-body perturbation theory 

In order to derive a microscopic approach to the effective interaction and/or operator within the framework of perturbation 
theory, we need to introduce various notations and definitions pertinent to the methods exposed. In this section we review 
how to calculate an effective operator within the framework of degenerate Rayleigh-Schrodinger (RS) perturbation theory 

mini 

It is common practice in perturbation theory to reduce the infinitely many degrees of freedom of the Hilbert space 
to those represented by a physically motivated subspace, the model space. In such truncations of the Hilbert space, the 
notions of a projection operator P onto the model space and its complement Q are introduced. The projection operators 
defining the model and excluded spaces are defined by 

D 

P = ^|<I>,)($,|, (1) 

i=l 

and 

oo 

Q= (2) 

i=D+l 

with D being the dimension of the model space, and PQ = 0, P^ = P, = Q and P + Q = I. The wave functions 1$^) 
are eigenfunctions of the unperturbed Hamiltonian Hq = T + U, where T is the kinetic energy and U and appropriately 
chosen one-body potential, that of the harmonic oscillator (h.o.) in most calculations. The full Hamiltonian is then 
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rewritten as H = Hq + Hi with Hi = V — U , V being e.g. the NN interaction. The eigenvalues and eigenfunctions of 
the full Hamiltonian are denoted by |\I/a) and Ea, 

(3) 

Rather than solving the full Schrodinger equation above, we define an effective Hamiltonian acting within the model space 
such that 

where — P I'J'q) is the projection of the full wave function onto the model space, the model space wave function. In 
RS perturbation theory, the effective interaction HcS can be written out order by order in the interaction Hi as 

PH^sP = PHiP + PHi—HiP + PHi—Hi — HiP + . . . , (5) 

e e e 

where terms of third and higher order also include the aforementioned folded diagrams. Further, e = uj — Hq, where to 
is the so-called starting energy, defined as the unperturbed energy of the interacting particles.. Similarly, the exact wave 
function can now be written in terms of the model space wave function as 

= + -Hi + ^Hi^Hi + . . . (6) 

e e e 

The wave operator is often expressed as 

= 1 + X, (7) 

where x is known as the correlation operator. The correlation operator generates the component of the wave function in 
the Q-space and must therefore contain at least one perturbation. Observing that PUP = P, we see that the correlation 
operator x has the properties 

PxP = 0, QflP - QxP - xP (8) 

Since — ft j^'f^) determines the wave operator only when it operates to the right on the model space, i.e., only the 
fiP part is defined, the term flQ never appears in the theory, and we could therefore add the conditions QxQ = and 
PxQ = to Eq. This leads to the following choice for x 

X = QxP (9) 

This has been the traditional choice in perturbation theory ^1 EI ■ 

The wave operator fl can then be ordered in terms of the number of interactions with the perturbation Hi 

n = i + n'-^^ + n''^^ + (lo) 

where fi*-"^ means that we have n Hi terms. Explicitely, the above equation reads 

where e are the unperturbed energies of the P-space and Q-space states defined by Hq. The greek letters refer to Q- 
space states, whereas a latin letter refers to model-space states. The second term in the above equation corresponds to 
fi^^^ while the third and fourth define fi^^^ Note that the fourth term diverges in case we have a degenerate or nearly 
degenerate model space. It is actually divergencies like these which are to be removed by the folded diagram procedure 
to be discussed in the next section. Terms like these arise due to the introduction of an energy independent perturbative 
expansion. Conventionally, the various contributions to the perturbative expansion are represented by Feynman-Goldstone 
diagrams. In Fig. ^we display the topologically distinct contributions to first order in the interaction V to the wave 
operator Eq. (|12|l . These diagrams all define the correlation operator x to first order in V . Diagrams with Hartree-Fock 
contributions and single-particle potential terms are not included. The possible renormalizations of these diagrams will 
be discussed in the next four sections. 

The linked-diagram theorem jl7l I41| can be used to obtain a perturbative expansion for the energy in terms of the 
perturbation V ov V — H — Hq where Hq represents the unperturbed part of the Hamiltonian. The expression for the 
energy E reads 



i? = ^(v|/o|i/[(c.-i/o)-^i/]'|*o)i, (12) 

fe=0 
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(e) (f) (g) 



Figure 1: The various vertices to first order in the interaction V which contribute to the wave operator CI = 1 + %. 
Hartree-Fock terms are not included. Possible hermitian conjugate diagrams are not shown. Indicated are also possible 
angular momentun coupling orders. 
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Figure 2: Antisymmctrizcd Goldstone diagrams through third order in perturbation theory included in the evaluation of 
the binding energy. The dashed lines represents the interaction. Particle and hole states are represented by upward and 
downward arrows, respectively. The first order diagram is omitted. All closed circles stand for a summation over hole 
states. 
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where vI/q is the uncorrelated Slater determinant for the ground state, u is the corresponding unperturbed energy and the 
subscript L stands for hnked diagrams only. In Fig. El we show all antisymmetrized Goldstone diagrams for closed-shell 
systems through third order in perturbation theory (we omit the first order diagram). All closed circles stand for a 
summation over hole states. The expression for e.g., diagram (4) AE^ in Fig. breads in an angular momentum coupled 
basis (with J being the total two-body angular momentum) 

Ai?4 = i {2J+l)mJ\V\{pq)J) ^ (13) 

ij < F 
pqrs > F 
J 

X {ipq)J\ V\{rs)J) [ _ ((rs) J| V \{ij)J) . 

£"2 I Sj Sj' Sg 

The symbol F represents the last hole, or Fermi surface. In the next section we will replace the bare interaction with the 
so-called G-matrix. 

We end this section with the equations for the diagrams in Fig. ^ representing x to first order in V. Moreover, in 
order to introduce the various channels needed to sum the parquet class of diagrams, we will find it convenient here 
to classify these channels in terms of angular momentum recouplings. Later on, we will also introduce the pertinent 
definitions of energy and momentum variables in the various channels. The nomenclature we will follow in our labelling 
is that of Blaizot and Ripka, see Ref. 0Tj chapter 15. All matrix elements in the definitions below are antisymmetrized 
and unnormalized. The first channel is the [12] channel, or the s-channel in field theory, and its angular momentum 
coupling order is depicted in Fig. 13 In this figure we do not distinguish between particles and holes, all single-particle 




(a) (b) (c) 

Figure 3: Coupling order for the [12] (a), [13] (b) and [14] (c) channels. 



labels 1,2,3,4 can represent either a hole or particle single-particle state, 
here. The matrix element V^I^^l in the [12] channel is then 



■12] 
1234J 



= ((12)J|y|(34)J) 



It is the coupling order which is of interest 



(14) 



meaning that the single-particle state 1(3) couples to the state 2(4) to yield a total angular momentum J. This channel 
is commonly denoted as the particle-particle (pp) or particle-particle-hole-hole (pphh) channel, meaning that when we 
will sum classes of diagrams to infinite order in this channel, the only intermediate states which are allowed are those of 
a pphh character, coupled to a final J in the above order. In the next section we will explicitely discuss ways to evaluate 
the equations for the [12] channel. This coupling order is also the standard way of representing effective interactions for 
shell-model calculations. If we now specialize to particles and holes (these matrix elements were shown in Fig. ^ we 
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obtain for the case with particle states only^, diagram (a), 

V2^^vl]^lj^{(j>q)J\V\{rs)J). (15) 
The corresponding one for holes only, diagran (b), is 

V2^^v}^^\,j^{{aP)J\V\{l5)J). (16) 

Thus, in the forthcoming discussion, we will always employ as our basic notation for a matrix element that of the [12] 
channel, meaning that matrix elements of the other two channels can always be rewritten in terms of those in [12] channel 
We see this immediately by looking at the expression for the matrix element in the [13] channel, the t-channel in field 
theory, see Fig.Ofb). Here the single-particle state 3(4) couples to the single-particle state 1(2)^. Through simple angular 
momentum algebra we have 

where the symbol with curly brackets represents a 6j-symbol and J' = \/2J' + 1. In a similar way we can also express 
the matrix element in the [14] channel, the u-channel in field theory, through 

vlllj = Y.(-r^''^'^'''j"{l I ^jW^sL- (18) 

It is also possible to have the inverse relations or to express e.g., the [14] channel through the [13] channel as 

w^L-E(-)'"^'^^^'">'{ j'Yi^^j- (19) 

The matrix elements defined in Eqs. (|14|I - H19|) and the inverse relations are the starting points for various resummation of 
diagrams. In the next section we will detail ways of solving equations in the [12] channel, whereas various approximations 
for the [13] channel and [14] channel such as the TDA and RPA and vertex and propagator renormalization schemes will 
be discussed in section O Finally, how to merge self-consistently all three channels will also be discussed in section |S| 

We end this section by giving the expressions in an angular momentum basis for the remaining diagrams of Fig. ^ 
The coupling order is indicated in the same figure. 

The 2plh vertex V2pih, diagram (c) in Fig.Q] is coupled following the prescription of the [13] channel and reads 

^2pi. = Vlfi, = X^(-)-+>+.^+-' J'^ { > I 'j, } Vlfi,. (20) 

The 2p2h ground-state correlation V2p2h, diagram (d), which will enter in the RPA summation discussed in section [S] is 
given by, the coupling order is that of the [13] channel, 

Vo - - \"(-V^+ip+J+J' / -^'^ \ T/[i2] (21) 

l'2p2h — ypqapj — ) \ jq j p J' J P1°'PJ" ^ ' 

The 2hlp vertex V2hip, diagram (e), still in the representation of the [13] channel, is defined as 

V■,^.^ -^l\]..+3v + J+-J' v'^ { h 3c. J \ y[12] ^22) 

v'2hip - v„^^pj - 2^1 ) \ in jp J' j O'PlPJ'- ^ > 

Note well that the vertices of Eqs. H2U|) - (|22|) and their respective hermitian conjugates can all be expressed in the [14] 
channel or [12] channel as well. However, it is important to note that the expressions in the various channels are different, 
and when solving the equations for the various channels, the renormalizations will be different. As an example, consider 
the two particle-hole vertices V^h of Fig. i.e., diagrams (f) and (g). Diagram (g) is just the exchange diagram of (f) 
when seen in the [12] channel. However, if (f) is coupled as in the [13] channel, recoupling this diagram to the [14] 
channel will not give two particle-hole two-body states coupled to a final J but rather a particle-particle two-body state 
and a hole-hole two-body state. But why bother at all about such petty details? The problem arises when we are to 
sum diagrams in the [13] channel and [14] channel. In the [12] channel we allow only particle-particle and hole-hole 

^Noto that we only include angular momentum factors, other factors coming from the diagram rules ;58^, like number of hole lines, number 
of closed loops, energy denominators etc. are omitted here. 

■^In a Goldstone-Feynman diagram in an angular momentum representation, the coupling direction will always be from incoming single- 
particle states to outgoing single-particle states. 
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intermediate states, whereas in the [13] channel and [14] channel we allow only particle-hole intermediate states, else 
we may risk to double-count various contributions. If we therefore recouple diagram (f) to the [14] representation, this 
contribution does not yield an intermediate particle-hole state in the [14] channel. Thus, diagram (f), whose expression is 

^Ph - VpfiaqJ - ) 1 i/3 jq J' j PP'^lJ" ^ ' 

yields a particle-hole contribution only in the [13] channel, whereas the exchange diagram (g), which reads 

Vph - V^Baaj - / ,\ ) \ jp jq J' j P01°'J'' ^ ' 

results in the corresponding particle-hole contribution in the [14] channel. In electron gas theory, the latter expression is 
often identified as the starting point for the self-screening of the exchange term. In the discussion of the TDA series in 
section 01 we will give the expressions for the screening corrections based on Eqs. H23|l and (|24|l . 

An important aspect to notice in connection with the latter equations and the discussions in section [3 is that 

y[14] _ _T^[13] f^r) 
^pPqaJ — ^ppaqj' l^^i 

i.e., just the exchange diagram, as it should be. This is however important to keep in mind, since we later on will sum 
explicitely sets of diagrams in the [13] channel and the [14] channel, implying thereby that we will obtain screening and 
vertex corrections for both direct and exchange diagrams. 



3 Diagrams in the [12] channel 

In order to write down the equation for the renormalized interaction r^^^l in the [12] channel we need first to present 
some further definitions. We will also assume that the reader has some familiarity with the theory of Green's function. In 
our presentation below we will borrow from the monograph of Blaizot and Ripka |41j , see also the recent review articles 
of Dickhoff and Barbieri and Miither and Polls _^52s- The vertex pl^^l is in lowest order identical with the interaction 
^[^■^1 and obeys also the same symmetry relations as V, i.e., 

p[12] _ p[12] _ _-p[12] _ -p[12] /r,px 
^ 1234J ~ ^ 2143J ~ ^ 2134J ~ ^ 1243J- y^^) 

We also need to define energy variables. Since we are going to replace the interaction V with the G-matrix, or certain 
approximations to it, defined below in all of our practical calculations, the momentum variables are already accounted 
for in G. The basis will be that of harmonic oscillator wave functions, and the labels 1234 will hence refer to oscillator 
quantum numbers, which in turn can be related to the momentum variables. The labels 1234, in addition to representing 
single-particle quantum numbers, define also the energy of the single-particle states. With a harmonic oscillator basis, the 
starting point for the single-particle energies £1,2,3,4 are the unperturbed oscillator energies. When iterating the equations 
for self-consistent single-particle energies can be introduced. The total energy in the [12] channel s is 

s = £1 4- £2 = £3 + £4. (27) 
The equation for the vertex r^^^l is, in a compact matrix notation, given by |41| 

r[i2i = F[i2] + ^[i2](^^)r[i2], (28) 

where g is the one-body Green's function representing the intermediate states. The diagrammatic expression for this 
equation is given in Fig. 01 The expression of Eq. (|28|l is known as the Feynman-Galitskii equation. This equation 
is normally solved iteratively. In the first iteration the irreducible vertex yl^^l is then often chosen as the bare NN 
interaction. This interaction is then typically assumed to be energy independent and we can drop the s dependence of 
I/I^^l . Moreover, the single-particle propagators are chosen as the unperturbed ones. The first iteration of Eq. (|28|l can 
then be rewritten as 

^1234 J (■5) — W234J + 2 5Z W256J^'^^'^5634j(^)- (29) 
56 

We have defined the unperturbed particle-particle and hole-hole propagators 
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Figure 4: (a) represents the two-body vertex T function while (b) represents the self-energy E. 



which results from the integration over the energy variable in the product of the two single-particle propagators in Eq. 
(I28|l . In our discussions we will not deal with dressed propagators, neither for the one-particle nor the two-particle 
propagators. For this and related topics such as the spreading of single-particle strength we refer to the work of Dickhoff 
and Barbieri |56j . In our discussion of the coupled cluster method however, such features can be extracted at the end of 
the calculations through various expectation values. These expectation values will however depend on the actual size of 
the model space used. In our case this applies to the number of harmonic oscillator shells. 

The factor 1/2 in Eq. (|29|l follows from one of the standard Goldstone-Feynnian diagram rules which state that 
a factor 1/2 should be associated with each pair of lines which starts at the same interaction vertex and ends at the same 
interaction vertex. This rule follows from the fact that we sum freely over the intermediate single-particle states 56. The 
reader should note that the intermediate states 56 can represent a two-particle state or a two-hole state. In Eq. H3()|l we 
have assumed unperturbed single-particle energies. In our iterations we will approximate the single-particle energies with 
their real part only. Thus, the two-particle propagator Q^^'^^ with renormalized single-particle energies has the same form 

[121 [121 

as the unperturbed one. The operators Qpp and Ql^^i ensure that the intermediate states are of two-particle or two-hole 
character. In order to obtain a self-consistent scheme, Eq. (I29|l has also to be accompanied with the equation for the 
single-particle propagators g given by Dyson's equation 



g = 90 + 90^9, 



(31) 



with go being the unperturbed single-particle propagator and E the self-energy. We will however defer a discussion of 
these quantities to sectional Here it will suffice to state that the self-energy is related to the vertex fI^^J as 



E^gT. 



(32) 



The similarity sign is meant to indicate that, although being formally correct, great care has to be exercised in order not 
to double-count contributions to the self-energy 001 ■ The set of equations for the vertex function and the self-energy is 
shown pictorially in Fig. 0] Assume now that we have performed the first iteration. The question which now arises is 
whether the obtained vertex F^^^l from the solution of Eq. H29|) should replace the bare vertex yl^^l in the next iteration. 
Before answering this question, let us give some examples of diagrams which can be generated from the first iteration. 
These contributions are given by e.g., diagrams (a)-(d) in Fig. Diagrams (a) and (b) are examples of contributions 
to second order in perturbation theory, while diagrams (c) and (d) are higher order terms. Diagrams (e) and (f) are 
higher-order core-polarization terms, which can e.g., be generated through the solution of the equations for the [13] and 
[14] channels discussed in the next section. If we were to include diagrams (a)-(d) in the definition of the bare vertex 
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Figure 5: Diagrams (a)-(d) give examples of diagrams which are summed up by the use of Eq. (|28|l . Diagrams (e) and 
(f) are examples of core-polarization terms which are not generated by the [12] channel. 



in our next iteration, we are prone to double-count, since such contributions are generated once again. Diagrams which 
contain particle- hole intermediate state are however not generated by the solution of Eq. (|29|l . We need therefore to 
define the vertex yl^^l used in every iteration to be the sum of diagrams which are irreducible in the [12] channel. With 
irreducible we will mean all diagrams which cannot be reduced to a piece containing the particle states 12 entering or 
leaving the same interaction vertex and another part containing the states 34 at the same interaction vertex by cutting 
two internal lines. Clearly, if we cut diagrams (a) and (b) we are just left with two bare interaction vertices. Similarly, 
cutting two lines of an intermediate state in diagrams (c) and (d) leaves us with two second-order terms of the type (a) 
and (b) and one bare interaction. Diagrams (e) and (f) are however examples of diagrams which are irreducible in the [12] 
channel. Diagram (e) is irreducible in the [13] channel, but not in the [14] channel. Similarly, diagram (g) is reducible in 
the [13] channel and irreducible in the [14] channel. This means that, unless we solve equations similar to Eq. (|29|l in the 
[13] channel and [14] channels as well, changes from further iterations of Eq. (|29|l will only come from the single-particle 
terms defined by Dyson's equation in Eq. (|^ . 

In the remaining part of this section, we will try to delineate ways of solving the above equations, and discuss 
possible approximations, their merits and faults. First of all, we will reduce the propagator of Eq. (|5n|l to only include 
particle-particle intermediate states. This will lead us to the familiar C?-matrix in nuclear many-body theory. Based on 
the G-matrix, we will construct effective interactions through perturbative summations. Applications of such effective 
interacions to selected nuclei will then be discussed. Thereafter, we will try to account for hole-hole contributions and 
self-consistent determinations of the single-particle energies through the solution of Dyson's equation. 



3.1 The G-matrix 

In nuclear structure and nuclear matter calculations one has to face the problem that any realistic NN interaction 
V exhibits a strong short-range repulsion, which in turn makes a perturbative treatment of the nuclear many-body 
problem prohibitive. If the interaction has a so-called hard core, the matrix elements of such an interaction (V'l V 
evaluated for an uncorrelated two-body wave function 'ip{r) diverge, since the uncorrelated wave function is different from 
zero also for relative distances r smaller than the hard-core radius. Similarly, even if one uses interactions with softer 
cores, the matrix elements of the interaction become very large at short distances. The above problem was however 
overcome by introducing the reaction matrix G^^^l (displayed by the summation of ladder type of diagrams in Fig. |S1 
with particle-particle intermediate states only), accounting thereby for short-range two-nucleon correlations. The G^^^l- 
matrix represents just a subset to the solution of the equations for the interaction pl^^l in the [12] channel, we have 
clearly neglected the possibility of having intermediate states which are of the hole-hole type. The matrix elements of the 
interaction yl^^l then become 

(V'l Gli^l 1^) = (^1 ^[121 [VIZ) (33) 
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where is now the correlated wave function containing selected correlation from the excluded space. By accounting 
for these correlations in the two-body wave functon ^P, the matrix elements of the interaction become finite, even for 
a hard-core interaction V. Moreover, as will be discussed below, compared with the uncorrelated wave function, the 
correlated wave function enhances the matrix elements of V at distances for which the interaction is attractive. The type 
of correlations which typically are included in the evaluation of the G^^^l -matrix are those of the two-particle type. If we 

[121 

label the operator Q in this case by Qpp , we can write the integral equation for the G- matrix as 

^[12] ^ yll2] ^ yll2] g[12] (34) 

s — Ho + iri 

implicitely assuming that limr; — > oo. The variable s represents normally the unperturbed energy of the incoming 
two-particle state. We will suppress irj in the following equations. Moreover, since one is often interested only in the 
fj[i2]_]jja,trix for negative starting energies, the G^^^l-matrix commonly used in studies of effective interactions has no 
divergencies. Note also, that compared with Eq. (|29() . we express the G-matrix in terms of operators. The explicit form, 
with e.g., the sum over imtermediate states is implicitely included here. We can also write 

. ^[i2l+v-[^2lg[i^2i J Q^^^^G^^^Ks). (35) 

[121 

The former equation applies if the Pauli operator Qpp commutes with the unperturbed Hamiltonian Hq, whereas the 

[121 

latter is needed if [Hq, Qpp ] 7^ 0. Similarly, the correlated wave function is given as 

I*) = m + -^a^'^ w , (36) 

or 

I*) = IV') + QIp' ^[i2!^^[i2i Qlp'g'"' 1^) ■ (37) 
s — i^pp -tioWpp 

In order to evaluate the G^^^l -matrix for finite nuclei, we define first a useful identity following Bethe, Brandow and 
Petschek jSH]- Suppose we have two different G-matrices^, defined by 

G^^Vi + Vi—Gu (38) 
ei 

and 

G2 = 1^2 + ^^2 — G2, (39) 

62 

where Qi/ei and Qilei are the propagators of either Eq. or Eq. (p^ . Gi and G2 are two different G-matrices having 
two different interactions and/or different propagators. We aim at an identity which will enable us to calculate G\ in 
terms of G2, or vice versa. Defining the wave operators 

r2i = l + ^Gi, (40) 
ei 



and 

we can rewrite the above G-matrices as 
and 

Using these relations, we rewrite G\ as 



f^2-l + — G2, (41) 



62 



Gi^li^i, (42) 
G2 = V2^i- (43) 



Gi= Gx-G\{^x-\-^G^ + {^\-\-G\^^Gx 

Gt+Gt f^-^)Gi+r!^Gi-G^17i, (44) 



61 62 

^For notational economy, wo drop the superscript [^^1. Furthermore, in the subsequent discussion in this subsection it is understood that 
all operators Q refer to particle-particle intermediate states only. The subscript pp is also dropped. 
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and using Eqs. (|42() and H43() we obtain the identity 



Gi = Gt + Gl -^9i]g, + nliv, - V2)ni. (45) 

The second term on the rhs. is called the propagator-correction term; it vanishes if Gi and G2 have the same propagators. 
The third term is often referred to as the potential-correction term, and it disappears if Gi and G2 have the same 
potentials. The reader may now ask what is the advantage of the above identity. If we assume that by some physical 
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Figure 6: Classification of harmonic oscillator single-particle orbits. 



reasoning we are able to calculate G2 and that the expression for G2 can be calculated easily, and further that G2 is a 
good approximation to the original G-matrix, then we can use the above identity to perform a perturbative calculation 
of Gl in terms of G2. 

Before we proceed in detailing the calculation of the G-matrix appropriate for finite nuclei, certain approximations 
need to be explained. 

As discussed above, the philosophy behind perturbation theory is to reduce the intractable full Hilbert space problem 
to one which can be solved within a physically motivated model space, defined by the operator P. The excluded degrees 
of freedom are represented by the projection operator Q. The definition of these operators is connected with the nuclear 
system and the perturbative expansions discussed in sectional Consider the evaluation of the effective interaction needed 
in calculations of the low- lying states of ^*0. From experimental data and theoretical calculations the belief is that several 
properties of this nucleus can be described by a model space consisting of a closed ^^O core (consisting of the filled Os- 
and Op-shells) and two valence neutrons in the IsOd-shell. In Fig. |H|we exhibit this division in terms of h.o. sp orbits. 
The active sp states in the IsOd-shell are then given by the 0^5/2, 0^3/2 and lsi/2 orbits, labels 4 — 6 in Fig. The 
remaining states enter the definition of Q. Once we have defined P and Q we proceed in constructing the G-matrix and 
the corresponding perturbative expansion in terms of the G-matrix. There are however several ways of choosing Q. A 
common procedure is to specify the boundaries of Q by three numbers, ni, n2 and 713, explained in Fig. For ^^O 
we would choose (rii = 3, n2 = 6,71.3 = 00). Our choice of P-space implies that the single-particle states outside the 
model space start from the lp0/-shell (numbers 7-10 in Fig. El, and orbits 1, 2 and 3 are hole states. Stated differently, 
this means that Q is constructed so as to prevent scattering into intermediate two-particle states with one particle in 
the Os- or Op-shells or both particles in the IsOd-shell. This definition of the Q-space influences the determination of 
the effective shell-model interaction. Consider the diagrams displayed in Fig. |S1 Diagram (a) of this figure is just the 
G-matrix and is allowed in the definition of the Q-box. With our choice (ni — 3,n2 = 6,^3 — 00), diagram (b) is not 
allowed since the intermediate state consists of passive particle states and is already included in the evaluation of the 
G-matrix. Similarly, diagram (c) is also not allowed whereas diagram (d) is allowed. Now an important subtlety arises. 
If we evaluate the G-matrix with the boundaries (ni — 3,n2 — 10, 713 = 00), and define the P-space of the effective 
interaction by including orbits 4 to 6 only, then diagrams (b) and (c) are allowed if 7 < pi,p2 < 10 In this way we allow 
for intermediate two-particle states as well with orbits outside the model-space of the effective interaction. The reader 
should notice the above differences, i.e. that the Q-space defining the G-matrix and HoB may differ. In order to calculate 
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Figure 7: Definition of the P (shaded area) and Q operators appropriate for the definition of the G-matrix and the 
effective interaction. Outside the shaded area hmited by the boundaries n\, n2 and ns P = and Q = 1. 
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Figure 8: Examples of diagrams which may define the effective valence space interaction. The wavy line is the G-matrix. 
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the G-matrix we will henceforth employ a so-called double-partitioned scheme. Let us be more specific and detail this 
double-partitioned procedure. We define first a reference G-matrix G in terms of plane wave intermediate states only, 
meaning that Hq is replaced by the kinetic energy operator T only while G has harmonic oscillator intermediate states 
(this is one possible choice for U). We divide the exclusion operator into two parts, one which represents the low- lying 
states Qi and one which accounts for high- lying states Qh, viz. 

Q = Qi + Qh = Qi + Q. 

If we consider ^^O as our pilot nucleus, we may define Qi to consist of the sp orbits of the lpO/-shell, orbits 7 — 10 in 
Fig. El described by h.o. states. Qh represents then the remaining orthogonalized intermediate states. Using the identity 
of Bethe, Brandow and Petschek [5^1 of Eq. H45|) we can first set up G as 



G^V + V- 



T 



G, 



and express G in terms of G 



as 



G = G + G 



Qi 



Ho, 



G, 



(46) 



(47) 



and we have assumed that G is hermitian and that [Qi,Hq] — 0. Thus, we first calculate a "reference" G-matrix (G in 
our case), and then insert this in the expression for the full G-matrix. The novelty here is that we are able to calculate G 
exactly through operator relations to be discussed below. In passing we note that G depends significantly on the choice 
of Hq, though the low-lying intermediate states are believed to be fairly well represented by h.o. states. Also, the authors 
of ref. |60| demonstrate that low-lying intermediate states are not so important in G-matrix calculations, being consistent 
with the short-range nature of the NN interaction. Since we let Qi to be defined by the orbits of the lpO/-shell, and 
the energy difference between two particles in the sd-shell and pf shell is of the order —14 MeV, we can treat G as a 
perturbation expansion in G. Eq. H47|l can then be written as 



G = G + G 



Qi 



Ho 



G + G 



Qi 



Ho 



G 



Ho 



G 



(48) 



The only intermediate states are those defined by the lpO/-shell. The second term on the rhs. is nothing but the second- 
order particle-particle ladder. The third term is then the third-order ladder diagram in terms of G. As shown by the 
authors of ref. |6U| , the inclusion of the second-order particle-particle diagram in the evaluation of the Q-box, represents 
a good approximation. The unsettled problem is however how to define the boundary between Qi and Qh- 

Now we will discuss how to compute G. One can solve the equation for the G-matrix for finite nuclei by employing a 
formally exact technique for handling Q discussed in e.g., Ref. j6U| . Using the matrix identity 



1 



1 



with A = s — T, to rewrite Eq. H46|l as'' 
where Gf is the free G-matrix defined as 



Q = ---P- 
QAQ A A PA-^P' A 



G = Gf + AG, 
1 



G, 



v + v- 



Gf 



s~T 

The term AG is a correction term defined entirely within the model space P and given by 



AG 



^v^p^±^pU 



A PA-^P A 

Employing the definition for the free G-matrix of Eq. (|51|1 . one can rewrite the latter equation as 

1 ^ 1 



AG 



-Gf-P^ 



P{e-^ + e-^GFe-^)P e 



~:P-Gf. 



(49) 

(50) 
(51) 

(52) 
(53) 



with e = s — T. We see then that the G-matrix for finite nuclei is expressed as the sum of two terms; the first term is 
the free G-matrix with no Pauli corrections included, while the second term accounts for medium modifications due to 
the Pauli principle. The second term can easily be obtained by some simple matrix operations involving the model-space 

■^We will omit the label G for the G-matrix for finite nuclei, however it is understood that the G-matrix for finite nuclei is calculated 
according to Eq. 1461 This means that we have to include the particle-particle ladder diagrams in the Q-box. 
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matrix P only. However, the second term is a function of the variable n^. The convergence in terms of 713 was discussed 
ad extenso in Ref. \W\ and we refer the reader to that work. The equation for the free matrix G_f is solved in momentum 
space in the relative and centre of mass system and thereafter transformed to the relevant expression in terms of harmonic 
ocillator single-particle wavefunctions in the laboratory system. This yields final matrix elements of the type 

{{ab)J\G\{cd)J) (54) 

where G is the given by the sum G ~ Gp + AG. The label a represents here all the single particle quantum numbers 

3.2 Modified Hilbert space for resummation of large sets of diagrams 

In Fig. 13 the numbers ni, n2 and 123 where used to define hole and particle states with respect to a given nucleus. We 
could also define a huge model space which does not reflect a particular core and thereby nucleus. One possible way 
of defining such a no-core model space is obtained by setting ni = n2 = n^, with 713 representing a large number, at 
least some eight-ten major oscillator shells. The single-particle states labeled by na represent then the last orbit of the 
model space P, following the numbering indicated in Fig. EJ This so-called no-core model space is indicated in Fig. 
and will be used in our definitions of model spaces for the resummations of parquet diagrams and many-body terms in 
coupled cluster theory. In Fig. |3 the two-body state \{pq)JTz) does not belong to the model space and is included in 
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Figure 9: Definition of the exclusion operator used to compute the G-matrix for large spaces. 

the computation of the G-matrix. Similarly, \{p'-f)JTz) and \{Sq)JTz) also enter the definition of Q whereas \ {6j)JTz) is 
not included in the computation of G. This means that correlations not defined in the G-matrix need to be computed by 
other non-perturbative resummations or many-body schemes. This is where the coupled-cluster scheme and the parquet 
approaches enter. 

With the G-matrix model space P of Fig.|5|we can now define an appropriate space for many-body perturbation theory, 
parquet diagrams or coupled-cluster calculations where correlations not included in the G-matrix are to be generated. 
This model space is defined in Fig. 1101 where the label Up represents the same single-particle orbit as 713 in Fig. 13 

The G-matrix computed according to Fig. does not reflect a specific nucleus and thereby single-particle orbits which 
define the uncorrelated Slater determinant. For a nucleus like '*He the Osi/2 orbit is fully occupied and defines thereby 
single-hole states. These are labeled by Ua in Fig.^| For ^^O the corresponding hole states are represented by the orbits 
Osi/2, 0p3/2 and Opi/2- With this caveat we can then generate correlations not included in the G-matrix and hopefully 
perform resummations of larger classes of diagrams. 

To demonstrate the dependence upon the size of the enlarged model-space defined by rip. We present here results from 
third-order in perturbation theory for the binding energy of ^^O as function of the size of the model space and the chosen 
oscillator energy hui. These results are shown in Fig. 1111 

There are several features to be noted. First of all one notices that the results seem to stabilize between seven 
and eight major shells. The fact that the energies seem to converge at this level of truncation is a welcome feature 
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Figure 10: Definition of partiele and liole states for coupled-cluster, parquet diagrams and perturbative many-body 
calculations in large spaces. The orbit represented by Ua stands for the last hole state whereas Up represents the last 
particle orbit included in the G-matrix model space. The hole states define the Fermi energy. 
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Figure 11: Binding energy per particle E/A from third-order perturbation theory for "'^^O as fimction of the number of 
major harmonic oscillator shells N and the oscillator energy hw. For = 8 we have the optimal value of E/A = —7.12 
MeV at fiw = 13.6 MeV. The experimental value is E/A = -7.98 MeV. 
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which can be exploited in the coupled-cluster and parquet diagram calculations. These calculations, see below, are much 
more challenging from a computational point of view since we in principle generate a much larger class of diagrams, 
coupled-cluster calculations. 

Secondly, although the minimum shifts a little as function of the oscillator energy as we increase the oscillator space, 
we notice that as the number of major shells is increased, the dependence of the binding energy upon the oscillator 
parameter weakens. A similar feature is seen in the coupled-cluster calculations below. For ^^0 the minimum for seven 
shells takes place at EjA = -7.16 MeV for huj = 12.9 MeV and for eight shells we have EjA = -7.12 MeV at huj = 13.6 
MeV. For six shells we obtain EjA — —1 A1 MeV at Tiuj = 12.6 MeV. The curvature for larger values of TvjJ decreases 
with increasing number of shells TV. At Tuxi = 18 MeV we have for ^^^O and TV = 8 that d{E/A)/duj = 0.22, for iV = 7 we 
obtain d{E/A)/du! = 0.28 and iV = 6 we have d{E/A)/du) = 0.35. The reader should also note that in the limit hu Q 
we have £' 0. 

We do not expect to reproduce the experimental binding energies. This calculation does not include the Coulomb 
interaction or realistic three-body interactions. The latter will be discussed in section 



3.3 Folded diagrams and the effective valence space interaction 

Here we discuss further classes of diagrams which can be included in the evaluation of effective interactions for the shell 
model. Here we will focus on the summations of so-called folded diagrams. 

One way of obtaining the wave operator 17 is through the generalized Bloch equation given by Lindgren and Morrison 

[ni 

Ha]P = QHiflP - xPHiQP, (55) 
which offers a suitable way of generating the RS perturbation expansion. Writing Eq. H55(l in terms of Sl^"-* we have 

[n'-'\Ho]P = QHiP, (56) 

[n'-^\Ho]P = QHiil'-^^P-fl'-^^PHiP, (57) 
and so forth, which can be generalized to 

[17("',i/o]P = QHifl^'^-^^P - J2 ^^"'^"'''PHin^""^^P. (58) 

m— 1 

The effective interaction to a given order can then be obtained from see |17) . Another possibility is obviously the 

coupled-cluster method discussed below. 

Here we will assume that we can start with a given approximation to Q, and through an iterative scheme generate 
higher order terms. Such schemes will in general differ from the order-by-order scheme of Eq. (|55|l . Two such iterative 
schemes were derived by Lee and Suzuki |HJ . We will focus on the folded diagram method of Kuo and co-workers • 

Having defined the wave operator fl = 1 + x (note that fl^^ = I — x) with x given by Eq. ((HJ we can obtain 

QHP - xHP + QHx - xHx = 0. (59) 

This is the basic equation to which a solution to x is to be sought. If we choose to work with a degenerate model space 
we define 

PHqP = sP, 

where s is the unperturbed model space eigenvalue (or starting energy) in the degenerate case, such that Eq. (|59|l reads 
in a slightly modified form {H = Hq + Hi) 

[s - QHoQ - QH,Q)x - QH,P - xPHiP - xPHiQx, 
which yields the following equation for x 

Observing that the P-space effective Hamiltonian is given as 

H,s = PHP + PHx = PHaP + V^x), 
with Vcsix) = PHiP + PHiQxP, Eq. ^ becomes 

^^T^^^^^^-T^^^^-W. (61) 
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Now we find it convenient to introduce the so-called Q-box, defined as 

Q{s) = PH,P + PH,Q L^QH.P. (62) 

s - QHQ 

The Q-box is made up of non- folded diagrams which are irreducible and valence linked. A diagram is said to be irreducible 
if between each pair of vertices there is at least one hole state or a particle state outside the model space. In a valence- 
linked diagram the interactions are linked (via fermion lines) to at least one valence line. Note that a valence-linked 
diagram can be either connected (consisting of a single piece) or disconnected. In the final expansion including folded 
diagrams as well, the disconnected diagrams are found to cancel out ^^l- This corresponds to the cancellation of unlinked 
diagrams of the Goldstone expansion. We illustrate these definitions by the diagrams shown in Fig. Diagram (a) is 
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Figure 12: Different types of valence-linked diagrams. Diagram (a) is irreducible and connected, (b) is reducible, while 
(c) is irreducible and disconnected. 



irreducible, valence linked and connected, while (b) is reducible since the intermediate particle states belong to the model 
space. Diagram (c) is irreducible, valence linked and disconnected. It is worth noting that general form of the Q-box is 
the same as that of the G-matrix, or the equations of the [12] channel or those of the [13] and [14] channels to be discussed 
in section^ In Ref. 16 , the Q-box was defined to be the sum all diagrams to third order in the G-matrix. 
Multiplying both sides of Eq. I|61(l with PHi and adding PHiP to both sides we get 



which gives 



PH^P + PHix = PHiP + PH,Q L^QH.P - PH, L^^V.six), 

s—QHQ s — QHQ 

Voffix) = - PHij^^^j^xVM- (63) 

There are several ways to solve Eq. (|63f) . The idea is to set up an iteration scheme where we determine Xn and thus 

Kff(Xn) from Xn-i and T4ff (x«-i). For the mere sake of simphcity we write V^^'^ = V;ff(Xn)- 
Let us write Eq. H63|) as 

The solution to this equation can be shown to be |61| 

ra — 1 

Observe also that the effective interaction is V^^^ is evaluated at a given model space energy s. If V^'j^'^ = V^^~^\ the 
iteration is said to converge. In the limiting case n — > oo, the solution V^^'' agrees with the formal solution of Brandow 
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[T^ and Des Cloizeaux |^ 



^ to! ds'^ 

m=0 



{V}^^y. (65) 



Alternatively, we can generate the contribution from n folds by the following way. In an n-folded Q-box there are of 
course n + 1 Q-boxes. The general expression for an n-folded Q-box is then 



TOi! ds"'i TO2! ds"2 '"TOnlrfs™" 



Q^Q I Q + Q I Q I Q^...^ ^ —-^P—-^P..—-^PQ, (66) 



mim2 . . .mn 



where we have the constraints 



mi + 7712 + . . . = 71, 

TOi > 1, 

1712, ms, ■ ■ ■mn>0, 



and 



TOfe < 'T- ^ + 1- 

The last restriction follows from the fact that there are only 7i — fc + 1 Q-boxes to the right of fc*'' Q-box. Thus, it can at 
most be differentiated 71 — fc + 1 times. We have inserted the model-space projection operator in the above expression, 
in order to emphasize that folded diagrams have as intermediate states between successive Q-boxes only model-space 
states. Therefore, the sum in Eq. (|66|l includes a sum over all model-space states with the same quantum numbers such 
as isospin and total angular momentum. It is understood that the Q-box and its derivatives are evaluated at the same 
starting energy, which should correspond to the unperturbed energy of the model-space state. It is then straightforward 
to recast Eq. Hfi6|l into the form of Eq. I|64|l . 

Note that although Q and its derivatives contain disconnected diagrams, such diagrams cancel exactly in each order 
|15| . thus yielding a fully connected expansion in Eq. I|64|) . However, in order to achieve this, disconnected diagrams have 
to be included in the definition of the Q-box. An example is given by diagram (c) Fig. 1121 Such a diagram will generate 

a contribution to the first fold ~ —Q J Q which cancels exactly diagram (c) when all time-ordered contributions to 
this diagram are accounted for, see Ref. |15| for more details. It is moreover important to note in connection with the 
above expansion, that a term like Fi = QiQ actually means PQiPQP since the Q-box is defined in the model space only. 

Here we have defined Qm — ;^4s^' ^'^^ ^^^^ structure, only so-called folded diagrams contain P-space intermediate 
states. 

The folded diagram expansion discussed above yields however a non-hermitian effective interaction. This happens 
even at the level of the G-matrix. A hermitian effective interaction has recently been derived by Suzuki and co-workers 
|lHlll8ll()Hj through the following steps^. To obtain a hermitian effective interaction, let us define a model-space eigenstate 
I 6a) with eigenvalue A as 

D 

|6a) = V^L^MV'o) (67) 



and the biorthogonal wave function 



such that 



D 

I^A)-E^i'^l^")' (68) 



Q = l 



{bx\h^,) ^ 5x^,. (69) 

The model-space eigenvalue problem can be written in terms of the above non-hermitian effective interaction unperturbed 
wave functions 

D 

b[^^ {iJa\Ho + V + VQx 1^7) = Exbi^'' ■ (70) 

7=1 



The exact wave function expressed in terms of the correlation operator is 

|*a) = (1 + x)|V'a). (71) 
The part x IV'a) can be expressed in terms of the time-development operator or using the time- independent formalism as 



The reader who wishes more details can consult Refs. IT^O . 
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where Q is the exclusion operator. Note that this equation is given in terms of the BriUouin-Wigner perturbation 
expansion, since we have the exact energy E\ in the denominator. 

Using the normahzation condition for the true wave function we obtain 

(^^IM'a) NxSx^ = (V^l (1 + x^x IV'a) , (73) 

where we have used the fact that (V'tI X IV'a) = 0. Recahing that the time-development operator is hermitian we have 
that x^X is also hermitian. We can then define an orthogonal basis d whose eigenvalue relation is 

E(^/3|X^X|V'a)rf^=M^4, (74) 

a 

with eigenvalues greater than 0. Using the definition in Eq. (I72|l . we note that the diagonal element of 

(V-aI xh IV'a) = (V'aI p^Qj^^—^j^Q^p ' (75) 

which is nothing but the derivative of the Q-box, with an additional minus sign. Thus, noting that if 7 7^ A 

(vj/^lM/;,) = = (V'^IVa) + ii^^lx^X I^a) , (76) 

we can write x^X i'^ operator form as 

X^X = - E F") (^"1 1^") (^"1 - E 1^") F/3> • (77) 

Using the new basis in Eq. I|74|l . we see that Eq. (|73|l allows us to define another orthogonal basis h 

^ vJ-^\ \/ui + 1 ^ 



p V ^ y t^a < p 

where we have used the orthogonality properties of the vectors involved. The vector h was used by the authors of Ref. 
|63 | to obtain a hermitian effective interaction as 

,,,,, ^/J€TT{^.\ K^-^ |#) + MTii^^l v:^-^ \^p) 



v/m^ + 1 + ^/4 + 1 



where (her) and (nher) stand for hermitian and non-hermitian respectively. This equation is rather simple to compute, 
since we can use the folded-diagram method to obtain the non-hermitian part. To obtain the total hermitian effective 
interaction, we have to add the Hq term. The above equation is manifestly hermitian. Other discussion of the hermiticity 
problem can be found in Refs. |17l I27| . The remaining question is how to evaluate the Q-box. Obviously, we are not in 
the position where we can evaluate it exactly, i.e., to include all possible many-body terms. Rather, we have to truncate 
somewhere. Several possible approaches exist, but all have in common that there is no clear way which tells us where 
to stop. However, as argued by the authors of Ref. [30], there is a minimal class of diagrams which need to be included 
in order to fulfil necessary conditions. This class of diagrams includes both diagrams which account for short-range 
correlations such as the G-matrix and long-range correlations such as those accounted for by various core-polarization 
terms. The importance of such diagrams has been extensively documented in the literature and examples can be found 
in Refs. 16, ,4.8^ In Ref. I16j we included all core-polarization contributions to third-order in the G-matrix, in addition 
to including other diagrams which account for short-range correlations as well. 

In our all results presented in section^jwe employ a third-order Q-box using the CD-Bonn interaction model described 
in Refs. 64, 65 . Moreover, the procedure for obtaining a hermitian procedure discussed in Eq. (|79() is employed in 
subsection 14.91 In our coupled-cluster calculations we employ the recent chiral model Idaho- A of Entem and Machleidt, 
see for example Ref. |66| . 



3.4 Center of mass corrections 

Momentum conservation requires that a many-body wave function must factorize as ^'(r) = (/)(R)^(ri.ci) where R is the 
center-of-mass (CoM) coordinate and rj-ci the relative coordinates. If we choose to expand our wave functions in the 
harmonic oscillator basis, then we are able to exactly separate the center-of-mass motion from the problem provided that 
we work in a model space that includes all nhfl excitations. In our coupled-cluster calculations we have a Q operator 
that allows for all possible two-particle interactions within a given set of oscillator shells. This means that we are nhVl 
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incomplete in a given calculation so that our method of separation of the center-of-mass motion becomes approximate. 
For example, for ^He in four major oscillator shells, we can excite all particles to n = 12fin excitations, but we can only 
excite one particle to n = 3h,^l excitations. Thus, care must be taken when correcting for center-of-mass contamination in 
our calculations. We have taken a variational approach based on the the work of Whitehead et ai, [HZI- The idea is to add 
/3coM-ffcoM to the Hamiltonian, but with /3coM remaining fairly small. This minimizes the effects of the center-of-mass 
contamination on low-lying state properties, and partially pushes unwanted states out of the spectrum. If we were to use 
a large /3coM , we would find spurious states entering into the calculated low- lying spectrum due to the incompleteness of 
our model space. 

The CoM Hamiltonian is then 

p2 1 3 

^^coM = 77TTT + ^mAn'^Ti^ - -nn , (80) 



2MA 



where P — a P* ^^'^ ^ ~ a ^i)/^- The term i?coM can be rewritten as a one-body harmonic potential, and 

a two-body term that depends on both the relative and center-of-mass coordinates of the two interacting particles. The 
matrix elements for the two body terms may be found in Ref. j68j . Operationally, we add -ffcoM to our Hamiltonian 



H' — H + f3coMH( 



CoM 



(81) 



where we choose /3coM so that the expectation value of i?coM is zero 69 . This insures that our center-of-mass contami- 
nation within the many-body wave function is minimized. We also find that this procedure yields reasonable spectra (in 
a space of four major oscillator shells) for ^He [TU]. This approach is used in our coupled-cluster calculations in sectional 
and in the shell-model Monte Carlo calculations of properties of IsOd— lpO/-shell nuclei, see section0]and Ref. j^. In 
our shell- model calculations we employ the standard procedure of Lawson, see Ref. j68| . 

There are at least two other widely used ways of dealing with the CoM corrections. The philosophy is 

• One starts with a translationally invariant Hamiltonian, kinetic energy plus e.g., two-body interaction. From that 
one the Harmonic Oscillator field is introduced. This is then split up in a pure CoM term and a two-body term 
which depends on the radial distance. 

• One starts with the the harmonic oscillator as the one-body piece plus a two-body interaction. The one-body piece 
can be rewritten in terms of a pure CoM term and a two-body term which now depends on the relative distance 
and momentum. 

In deriving the CoM corrections, the following expressions are helpful. The CoM momentum is 

A 



and we have that 



meaning that 



i<j 



9m 



2m 2mA 



The last expression is explicitely translationally invariant. 
In a similar fashion we can define the CoM coordinate 



2mA 



(82) 



(83) 



(84) 



1 ^ 



(85) 



which yields 



^ 1 



i=l 



A 



(86) 



If we then introduce the harmonic oscillator one-body Hamiltonian 



i=l 



2m 2 



(87) 



21 



with Q the osciUator frequency, we can rewrite the latter as 



2mA 



2mA ^ 



or just 



Hho = Hco 



M 



2mA 



Pj) 



with 



2A 



i<j 



(89) 



Hco 



P2 mAn'^R'^ 



M 



2mA 2 

In shell model studies the translationally invariant one- and two-body Hamiltonian reads for an A-nucleon system. 



(90) 



H 



^ -<2 

^ 2m 

.i=i 



2mA 



A 



(91) 



where m is the nucleon mass and Yij the nucleon-nucleon interaction, is modified by including the harmonic oscillator 
potential through the following 



1=1 



2^ 



2A 



0. 



We can rewrite the Hamiltonian of Eq. H91|l as 

1 



A 

E 

1=1 



2m 



A 

E 

i<3 



mQ^ 
~2A 



in 



CoM- 



(92) 



(93) 



Shell-model calculations are carried out in a model space defined by a projector P. The complementary space to the 
model space is defined by the projector Q = 1 — P. Consequently, for the P-space part of the shell-model Hamiltonian 
we get 



= E^ 



pI_ 

2m 



1 



P 



E^ 

i<j 



2A 



P - PHcomP. 



(94) 



The effective interaction appearing in Eq. (|94|l is in general an A-body interaction and if it is determined without 
any approximations, the model-space Hamiltonian provides an identical description of a subset of states as the full-space 
Hamiltonian (jHSl- The intrinsic properties of the many-body system still do not depend on fi. From among the eigenstates 
of the Hamiltonian JMJ, it is necessary to choose only those corresponding to the same CoM energy. This can be achieved 
by projecting the CoM eigenstates with energies greater than ^Tifl upwards in the energy spectrum. 

The effective interaction should be determined from (|93|l . Calculation of the exact A-body effective interaction 
is, however, as difficult as finding the full space solution. Usually, the effective interaction is approximated by a two-body 
effective interaction determined from a two-nucleon problem. The relevant two-nucleon Hamiltonian obtained from (|93|) 
is then 



ttVI 
^02 



p1 



p1 



2m +-,mn^iri+ri) + Vin-r,) 

With this Hamiltonian we can then compute a starting-energy-dependent effective interaction or G-matrix correspond- 
ing to a two- nucleon model space defined by the projector P2. This G-matrix reads 



-{ri-r2y 



(95) 



1 



(96) 



where Q2 — ^ — P2 and Vp is the interaction given by the last two terms on the righthand side of Eq. (|95|l . Using a 
similarity transformation introduced by Suzuki and Okamoto Navratil et al. used the above Hamiltonian to define 
an effective two-body interaction for no-core calculations, see for example Ref. [T^. The CoM corrections can easily be 
implemented within such a scheme. For a standard G-matrix this is more difficult since such calculations are typically 
done in momentum space. We have however applied the CoM correction to the bare potential in a perturbative way in 
Ref. d]. 

If we however start with a harmonic oscillator basis we have 



A r 

^" = E 

i=l 



P- + i-mn'fj 

2m 2 



A 

Kj 



(97) 
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Table 1: Shell- model dimensionalities for selected nuclei around A 100. 



bystem 


Dimension 


System 


Dimension 


^'J^Sn 


w 1.5-10^ 


^^^Sn 


w 6.2 • 10« 




w 3.2- IQS 


ii^Sn 


w 1.6- 10^ 




« 6.5-103 


"2Sb 


w 1.1 • 10*^ 


lo^Sb 


w 3.2- 10^ 


ii^Sb 


w 1.9- 10^ 



If we then use Eq. H88fl . we obtain the following expression for our Hamiltonian 

^ = ^ + — ^ + ^ - p.)' + ^ E(-' - + E ' (98) 



2mA 2 2mA -^'^ 2A 

or 



= i?coM + j:iP. - ft)^ + ^ E(-^» - -1)' + E (99) 



2mA 2A 

i<j i<j i<j 



which differs from Eq. (|94|l by the sign of the CoM contribution and the two-body part which now reads 

mO^ 1 

V- = V{n ~ r2) - -j^iri r^f ^(pl - Vi)\ (100) 



which should be contrasted to 



A thorough discussion of various approaches to the CoM problem within the context of the coupled cluster method 
and summation of parquet diagrams will be presented in future works |72| . 



4 Shell-model studies of selected nuclei 

The shell model problem requires the solution of a real symmetric n x n matrix eigenvalue equation 

H\veck) = Ek\veck) , (101) 

with k = I,. . .,K. The eigenvalues Ek are understood to be numbered in increasing order. In a typical shell model 
problem we are interested in only the lowest eigenstates of Eq. p01|) . so K may be of the order of 10 to 50. The total 
dimension n of the eigenvalue matrix H is large, for the Sn isotopes of interest up to n « 2 x 10^. Consider for example 
the case of a shell model calculation outside the Z = 50 = 50 core with the single-particle orbits 1^5/2, O57/2, 1^3/21 
2s]^/2 O/111/2 defining the shell-model space. Tabled shows examples of dimensionalities encountered for different 
nuclei. Such large matrix problems are increasingly used in science and enginering and practical numerical algorithms 
for determining their properties are continually being developed. These methods are closely related to the development 
of modern computer technology, both hardware and software. What was impossible to solve a few years ago may now 
be within reach and we should also be prepared for an increased future development. To indicate the present possibility, 
in a work by J. Olsen et al. |73] in a quantum chemistry configuration interaction calculation, Eq. (|101|l was solved in 
a basis with n — 10^. In nuclear physics there are presently several groups dedicated in the development of shell-model 
codes which can deal with systems with similar large dimensionalities. Notably are the lowa-Livermore- Tucson group 
[731, the Michigan group with e.g., the OXBASH code |31|731|7n|, the Oak Ridge group |77], the Oslo group |7H1 and 
the Strasbourg group [7^], just to mention a few. All these codes are nowadays capable of dealing with nuclear systems 
with dimensionalities of n ^ 10^ basic states. Extensive full /p-shell model and partial sd — fp calculations have been 
carried out by the Strasbourg group |80l 1811 IH^ I83| . Similarly, the lowa-Livermore- Tucson group has dealt with no-core 
shell- model calculations of light nuclei Ej El El , with a very good agreement with the ab-initio Green's function Monte 
Carlo calculations of the Argonne group |10l 1111 IT^ . 

Even larger shell- model dimensionalities can be dealt with e.g., Monte Carlo shell- model approaches developed by 
Otsuka and collaborators [HI I84L 155] and the shell-model Monte Carlo approach of Koonin et al. 7 . 

Below we will focus on results for medium heavy nuclei in the mass regions from A = 100 to A = 140 based on the 
Oslo shell-model code [TS] . 

Different computational approaches to solve Ea. (|101|l can be distinguished based on the size of n. For n small, 
i.e. 10^ < n < 10'^ and with the number of matrix elements of H less than 10^ such problems can be accomodated 
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within the direct access memory of a modern work station and can be diagonahzed by standard matrix routines. In a 
second domain with n > 10"^ but small enough that H has no more than lO'^ elements. This will require ~ 1.5 Gbyte 
of storage. Then all matrix elements may be stored in memory or alternatively on a standard disk. In these cases the 
complete diagonalization of H will not be of physical interest and efficient iteration procedures have been developed to 
find the lowest energy eigenvalues and eigenvectors. 

Based on the present computer methods we have developed a code which is under continuous improvement to solve 
the eigenvalue problem given in Eq. H101|l . The basic requirement is to be able to handle problems with n > 10^. In the 
following we discuss some of the important elements which enter the algorithm. 

We separate the discussion into three parts: 

• The m-scheme representation of the basic states. 

• The Lanczos iteration algorithm. 

• The Davidson-Liu iteration technique. 



4.1 The m-scheme representation. 

We write the eigenstates in Eq. (|101|l as linear combinations of Slater determinants. In a second quantization represen- 
tation a Slater determinant (SD) is given by 

\SD4N))= n (102) 

and the complete set is generated by distributing the N particles in all possible ways throughout the basic one-particle 
states constituting the P-space. This is a very efficient representation. A single 15*1?) requires only one computer word 
(32 or 64 bits) and in memory a \SD) with N particles is given by 

\SD) — > ( pOllllOlOlO- ■; ), (103) 

Nl's 

where each and 1 corresponds to an m-orbit in the valence P-space. Occupied orbits have a 1 and empty orbits a 0. 
Furthermore, all important calculations can be handled in Boolean algebra which is very efficient on modern computers. 
The action of operators of the form al^ap or al^aj^a^as acting on an \SD) is easy to perform. 

The m-scheme allows also for a straightforward definition of many-body operators such as one-, two- and three-particle 
operators 

alap, (104) 
ai,al^ap,a,32, (105) 

ail ^2 43^/31^/32 a/33 > (106) 
respectively, or generalized seniority operators. The seniority operators can be very useful in preparing a starting vector 
for the Lanczos iteration process. This option is fully implemented in our codes. 
The generalized seniority operators '86' can then be written as 



for seniority zero. 



J ' rn>0 



dIm= E (l + '5w')"'/'/3„-'(jmj'm'|JM)a]„a],„, (108) 

for seniority two. The coefficients Cj and (3jji can be obtained from the a chosen two-particle system such as the ^^°Sn 
ground state and the excited states, respectively. We can also define a seniority four operator 

G(ni,ji,n2,j2;J,M) = {<,,i<,,2 

= E 9i''.u,<^lAA. (109) 



and a seniority six operator 

I{ni,ji, {n2 , j2 , ns , js ) j23 \J,M) = { An ji ^(^2 , 32 , ^3 , h ; ^23 ) } jm=o 



E 9i^^..>^AA2°'lAAAe (110) 



1^1 ...1/6 



Finally, our shell- model code allows also for the inclusion of effective and real three-body interactions. Results from 
such calculations will be discussed in section [7| 
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4.2 The Lanczos iteration process. 

At present our basic approach to finding solutions to Eq. (|1(J1|I is the Lanczos algorithm. This method was already applied 
to nuclear physics problems by Whitehead et al. in 1977. In a review article ,6 7, they describe the technique in detail. 
For the present discussion we outline the basic elements of the method. 

1. We choose an initial Lanczos vector \lanco) as the zeroth order approximation to the lowest eigenvector in Eq. (|101|l . 
Our experience is that any reasonable choice is acceptable as long as the vector does not have special properties 
such as good angular momentum. That would usually terminate the iteration process at too early a stage. 

2. The next step involves generating a new vector through the process \newp+i >= H\lancp >. Throughout this 
process we construct the energy matrix elements of H in this Lanczos basis. First, the diagonal matrix elements of 
H are then obtained by 



{lancpl H \lancp) — {lancp\ neWp+i) , 
3. The new vector \newp+i) is then orthogonalized to all previously calculated Lanczos vectors 



(111) 



neWp_^_i^ — |neWp_|_i) — \lanCp) ■ {lanCp\ neWp+i) — \lanCq) ■ {lanCq\ neWp+i) , 

9=0 



and finally normalized 



(112) 



\lancp+i) = 

to produce a new Lanczos vector. 
4. The ofF-diagonal matrix elements of H are calculated by 

{lancp+i\ H \lancp) — ( 



(113) 



new. 



p+i 



(114) 



and all others are zero. 
5. After n iterations we have an energy matrix of the form 



Hq^O 


Hq.i 








Hq.i 












H2,l 













Hp- 








Hpp 


-1 Hp 



(115) 



as the p'th approximation to the eigenvalue problem in Eq. H101|l . The number p is a reasonably small number 
and we can diagonalize the matrix by standard methods to obtain eigenvalues and eigenvectors which are linear 
combinations of the Lanczos vectors. 

6. This process is repeated until a suitable convergence criterium has been reached. 

In this method each Lanczos vector is a linear combination of the basic \SD) with dimension n. For n w 10^ — 10®, 
as in our case of interest. Here is one of the important difficulties associated with the Lanczos method. Large disk 
storage is needed when the number of Lanczos vector exceeds « 100. Another difficulty is found in the calculation of 
\neWp+i >= H\lanCp > when n > 10^. 

One important objection found in the computer literature |87| to the Lanczos method is its slow convergence. This is 
also our experience so far and means that a large number of Lanczos vectors have to be calculated and stored in order to 
obtain convergence. One improvement which we have implemented is to terminate the Lanczos process earlier, diagonalize 
the energy matrix and choose some of the lowest thus obtained eigenvectors as a starting point for a new Lanczos iteration. 
This modifies the energy matrix in Eq. H115|) slightly, since the matrix will not be tri-diagonal any more. It reduces the 
disk storage requirement, but the convergence problem is left. A possible way out here is the Davidson-Liu method. 
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4.3 The Davidson-Liu method. 



We outline the basic elements of this technique and refer to the literature [87' for details. The first important improvement 
to the Lanczos process given in Eqs. ()111I115|) is to start with several orthogonalized and normalized initial vectors and 
the second difference is the way new additional vectors are chosen. The method can be viewed as an improvement of the 
Lanczos technique and can be described in the following steps: 



1. Choose a (small) set of start vectors 



k = l. 



2. K new vectors are generated through the process 

3. The matrix elements of H are again obtained by 



H 



new, 



(0) 



^ V fc — 1, 



(116) 
K and K 



and the matrix H is no longer tri-diagonal. H is diagonalized within the set of states 
eigenvalues and corresponding eigenstates are obtained as the zero'th approximation to the lowest true eigenstates 
of Eq. ifTUT)) in the form 



K 
9=1 



(117) 



The vectors 



now constitute the basis for a new iteration. 



4. New correction vectors are calculated through the steps 



lip 



\S,) - {H,,a, - s(-^)-^ \rp) . (118) 
Then new additional vectors are obtained by orthogonalizing to all previous vectors, and finally normalized 
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(«) 



(119) 



witk k' — I, . . . , K' < K . Thus up to K new vectors may be generated through one iteration. Then the diagonal- 
ization process in Eqs. (|116H117|) is repeated and a new iteration to the true eigenvectors is obtained. 

5. Again this process may be repeated until some convergence critcrium has been reached. 



4.4 Break-Up of the Doubly-Magic Sn Core 

Doubly-magic nuclei and their immediate neighbors are of great interest as they provide direct information on the basic 
shell structure that is ultimately responsible for most nuclear properties. Lying at the proton drip line and being the 
heaviest particle-stable, self-conjugate nucleus, ^'"^Sn is particularly relevant in this context. An important property of this 
nucleus is the degree of rigidity of its spherical shape which is reflected in the excitation energy of the lowest 2+ state and 
in the associated B(E2;2+ ^0^) transition rate. The main component of the wave function of this level in a microscopic 
description is presumably an isoscalar mixture of proton and neutron 2^5/2159^2 excitations across the N—Z— 50 shell 
gaps. This state is at present not known experimentally and its observation may well require the availability of intense 
exotic beams. Some guidance about its excitation energy can perhaps come from other doubly-magic nuclei. In the 
N=Z =28 doubly-magic nucleus, the first 2+ state is located rather low, at 2.7 MeV. In contrast, the 2+ levels in ^"^^Sn 
and ^°*Pb are much higher in excitation energy, 4.0 and 4.1 MeV, respectively, and in the latter nucleus the size of the 
shell gaps can also be appreciated from the fact that this state is not even the lowest excitation, but instead lies above a 
3^ (octupole) vibrational state. 

In order to estimate the position of the 2+ state in ^"^Sn, both proton and neutron shell gaps have to be known. 
The energy splittings of the relevant single particle orbitals in the other heavy doubly magic nuclei are comparable. The 
neutron 2^3/2 and I/7/2 orbitals are 6.4 MeV apart in ^^Ni, while the splittings between the 2/7/2 and I/111/2 levels in 
^^^Sn and 2^9/2 and 1^13/2 states in ^°®Pb are 4.9 and 5.1 MeV, respectively. Here, the splitting between the 2^5/2 and 
1^9/2 neutron orbits will be shown to be of the order of 6 MeV as in the ^^Ni case. However, a sizable proton- neutron 
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Q-space: Z > 50 






Q-space: iV > 82 H| 







/in/2 3.5 MeV 

97/2 2.63 MeV 

0.9 MeV 2.23 MeV 

P-space for particles P-space for particles 

1-26 MeV 



Pi/2 0.00 MeV 0.00 McV 




Figure 13: Possible shell-model space for nuclei around A ^ 100 using ®*Sr as closed-shell core with the neutron orbitals 
2si/2, 1(^5/2, 1^3/21 0^7/2 and O/111/2 and proton single-particle orbitals Ogg/2 and lpi/2 defining the model space. Their 
respective single-particle energies are displayed as well. 



interaction could, as in ^^Ni, decrease the excitation energy and increase the transition rate for the lowest 2+ state in 
^''"Sn. Such an interaction is expected to be especially strong in N—Z nuclei, where protons and neutrons occupy the same 
single-particle orbitals, which results in a large spatial overlap of their wave functions. Nuclei near ^'^"Sn were studied 
using the ^^Ni-f '"^"Cr reaction. A ^*Ni beam of 225 MeV was provided by the ATLAS superconducting linear accelerator 
at Argonne National Laboratory. The effective two-body interactions were applied in a shell-model space spanning 2si/2j 
1^5/2, 11^3/2: O57/2 and O/111/2 neutron, and OgQ/2 and lpi/2 proton single-particle orbitals 85". The model space is shown 
in Fig. El In this model space, ^^Cd has one neutron and 10 protons outside the core. The results of the calculation, 
denoted as SMH, are compared with the experimental levels in Fig. The calculation favors a J^=5/2+ assignment 
for the ground state, in agreement with the systematics of odd-A, A^=51 isotones. The wave functions of the states with 
J^=5/2+, 9/2+, 13/2+, 17/2+, 19/2'^ and 21/2^ all have main configurations where the valence neutron is in the ^5/2 
orbit, while the 7/2+, 11/2+, 15/2+, 19/2^, 21/2^ and 23/2+ levels are associated with the predominant occupation 
of the gj/2 orbit. For all these states, the two proton holes always remain assigned to the gg/2 orbit. The lowest 7/2+ 
states in N=51 isotones from ^^Zr to ^^Pd are well reproduced by using the experimental 177/2 single-particle energy from 
zirconium [SS], but in ^^Cd this state was calculated 250 keV too low. Therefore, the single-particle energy of the neutron 
(/7/2 orbit, relative to the ^^Sr, core was increased from 2.63 to 2.89 MeV in order to reproduce the experimental excitation 
energy of the 7/2+ state. This yielded excellent agreement between calculations and experiment up to the ^^^=23/2+ 
state, see for example Fig. the highest spin that can be generated in this model space for positive-parity states in 
^^Cd. The description of higher-spin states requires the excitation of one or more gg/2 particles across the N, Z=50 shell 
gaps, an excitation similar to that responsible for the first 2+ and higher-lying states in -'^°°Sn. 

To study the high-spin states, another shell-model calculation was performed using this time ''^Sr as a closed shell 
core. The results are denoted as SMG in Fig. The same model space as in the previous calculation was employed 
except for the /111/2 neutron orbit which was replaced by the 179/2 one. (Since this /111/2 state lies at a relatively high 
excitation energy, its contribution to the positive-parity states is small.) Due to limitations in computing time, only up 
to two neutrons were allowed to leave the 39/2 orbit. (Note that the (79/2 orbit lies below the iV=50 shell gap, while all 
other neutron orbitals in this model space lie above this gap.) Opening of the iV=50 shell only is justified since a neutron 
hole in the gg/2 orbit, together with an existing proton hole, produces a very attractive 9+ proton-neutron coupling. 
Similarly, the excited neutron, together with the one already occupying the same orbit, produces a 0+ neutron pair with 
a strong attractive interaction. Opening the Z=50 shell would not result in such a strong attraction. Single-particle 
energies with respect to the ^*Sr core are not known and were kept the same as in Ref. [HE! for the *^Sr core. The energy 
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Figure 14: Experimental (EXP) and calculated (SMH, SMG) level schemes for ^^Cd. The widths of the arrows are 
proportional to the intensities of the 7 rays observed in the experiment. As discussed in the text, the J = 25/2+ and 
higher-lying states require excitations across the ^^^Sn core. 
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of the (79/2 orbit was placed 5.0 MeV below that of the one. The wave functions of the states below the 25/2+ 
level are very similar to those obtained in the SMH calculation (which assumed that the /2 neutron orbit is completely 
filled), except for the two 19/2+ levels, where the occupation numbers are reversed for the and 37/2 orbits. The 
25/2+ level and the higher- lying states have 9 neutrons in the 39/2 orbit with the remaining neutron pair almost evenly 
distributed over the and 37/2 orbits. These states represent the break-up of the doubly-magic ^''^Sn core. Thus, 
the excitation energy of the states with J'^>25/2+ is sensitive to the position of the 179/2 orbit. By fitting this single 
particle energy to 5.0 MeV we, therefore, indirectly deduced the size of the iV=50 shell gap to be 6.5 MeV as defined 
by 2BE(^"°Sn)-BE(^^Sn)-BE(^°^Sn), where BE stands for binding energy. This value agrees well with earlier predictions 
from Hartree-Fock calculations by Leander et al. [SS], as well as a single-particle energy parametrization by Duflo and 
Zuker |^ and an extrapolation by Grawe et al. [51] . To investigate the effect of the truncation of the model space on 
the deduced size of the A^=50 shell gap we calculated the excitation energy of the J'^=25/2+ state in ^^Cd using two 
different model spaces. In the first one, we allowed up to four neutrons to leave the ^9/2 orbit and the second truncation 
scheme restricted the valence particles to only two in each active orbit. Different truncations of the model space required 
adjustments of the (79/2 single-particle energy to reproduce the experimental excitation energy of the 25/2+ level, but 
when this was achieved the size of the iV=50 shell gap remained within 0.5 MeV of the above value. This illustrates a 
relative insensitivity of our inferred result to the truncation procedure used in the shell model calculation. 

Calculations using the same single particle energies and matrix elements were also performed for ^"^In, which has one 
proton hole and two neutrons outside the N=Z='bQ doubly-closed shells. Again, the results, labeled as SMH and SMG, 
are in very good agreement with the experimental level scheme in Fig. 1151 In these calculations, the proton hole remains 
in the 59/2 orbital and a J'^=9/2+ assignment follows for the ground state. This is in agreement with the systematics 
of odd-A In isotopes. The valence neutron pair is distributed over the (^5/2 and 57/2 orbits for the ground state, while 
it occupies predominantly the ^5/2 state for the 13/2+ level. One neutron is promoted to the 57/2 and /111/2 orbitals 
to generate the 17/2+, 19/2+, 21/2+ positive-parity, and 23/2^, 25/2~ negative-parity states, respectively. The large 
separation between the 19/2+ and 21/2+ states, which are associated with the same dominant configuration, is attributed 
to the strong repulsion between the aligned 59/2 proton hole and the 57/2 neutron. Note that the 7r(79/2J^/iii/2 matrix 
elements used to calculate the negative-parity states in ^°^In described in Ref. 1^2] were adopted here as well. The 
maximum spin that can be reached with the tt (79/ 2 1^^5/2 57/ 2 configuration is 21/2. Excitation of one or more particles 
across the iV, Z=50 shell gaps is again needed to account for the positive-parity states of higher spin. Thus, the 23/2+ 
level is the lowest-lying core-excited state in ^"^In. The calculated energy gap of 2941 keV between the 23/2+ and 19/2+ 
states is in excellent agreement with the experimental value of 2924 keV. Experimentally known lowest lying core excited 
states were also calculated in the nuclei ^*Ag [HSI and ^^Pd 94 using the SMG model space. Their excitation energies 
were reproduced to better than 100 keV. It is more difficult to identify core excited states in nuclei near ^°Zr, because in 
those nuclei high-spin states can be easily reached by promoting a pair of protons from the low-spin Pii2 orbit into the 
empty ^9/2 orbit. 

The interactions used in these SMH and SMG calculations were also applied to a shell model calculation of the lowest- 
lying levels in ^°^Sn. Both SMH and SMG calculations favor 7^^=5/2+ quantum numbers for the ground state, while 
the 7/2+ state lies only ~100 keV above. This is in excellent agreement with the extrapolated energy deduced from the 
systematics of odd-A Sn isotopes down to ^°^Sn ^5) . 

It is also worth pointing out that the level schemes of ^^Cd and ^°^In closely mirror those of their analogs in the ^^Ni 
region; i.e., ^^Fe jHEj and ^^Co [HZI- In particular, the lowest-lying core-excited states have almost identical excitation 
energies. This observation, together with the deduced iV=50 shell gap, again points to a large similarity between the two 
heaviest self-conjugated doubly-magic nuclei and one may conclude that the lowest 2+ states in ^^Ni and ^^''Sn would 
have closely related excitation energies if the sizes of the Z=28 and Z=50 shell gaps were also similar. 

4.5 ^^^Sn Core Excitations in ^^"^In 

In the above model space ^"^In has three valence neutrons and 11 valence protons, which is equivalent to one proton 
hole in the doubly magic ^""^Sn core. In all calculated levels shown in Fig. 1161 this proton hole remains in the (79/2 orbit. 
Therefore, the only proton contribution to the calculated level scheme of ^''^In is through proton-neutron pairing. The 
calculation favors J'^— 6+ for the ground state. The three valence neutrons are mainly in the d^/2 orbit in this state, 
while about 40% of the wave function amplitude comes from the contribution of the 57/2 orbit. All other orbits have an 
insignificant contribution to the wave function of the ground state. Surprisingly, all levels up to the J'^=10j level have a 
very similar wave function configuration in the calculation. The 10^ level is the lowest observed level for which the ^5/2 
and (77/2 orbits switch their occupation numbers. The larger wave function difference results in lower mixing between the 
10^ and 103" levels. This in turn leads to two close lying 10+ states, correctly predicted by the calculation. However, all 
other non-yrast states are calculated too high and we expect that the (77/2 orbit has a larger than calculated contribution 
to their wave functions. This is especially true for the 11 J state that is calculated more than 1 MeV too high. We may 
conclude that the interaction used in the calculation gives a too strong attraction between the gg /2 protons and the 57 /2 
neutrons. 
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Figure 15: Experimental (EXP) and calculated (SMH, SMG) level schemes for ^°^In. The widths of the arrows are 
proportional to the intensities of the 7 rays observed in the experiment. The (J = 23/2+) state requires excitations 
across the ^"^Sn core. 
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Figure 16: Comparison of experimental and calculated levels in ^"^In. 



The very nice one to one correspondence between experimental and calculated positive parity states ends with the 
13^ state, which is the highest spin that can be reached by coupling one proton hole in the 179/2 orbit with three neutrons 
in c?5/2 and 37/2 orbits. The wave function of this state is therefore 7r((79/2)~^i^(i5/2(57/2)^- 

The experimental positive parity levels above the 3858 keV level do not have calculated counterparts. The only 
possibility to reach such spins and parities within the model space used in the calculation is to promote a neutron pair 
to the /111/2 orbit. However, even after changing the 7r59/2i^^ii/2 effective interaction to reproduce the negative parity 
states the J^=12>2 to l?"*" levels were calculated much higher in energy than the experimental ones. These states must, 
therefore, be due to excitations across the doubly closed N — Z — shell. Most likely they are due to the excitation of 
a neutron from the (79/2 orbit just below TV = 50 shell gap to the ^5/2 orbit just above the gap, since this produces the 
very attractive {'^9^J^2^gQ/^2T~^ coupling combined with four neutron particle states like the J'^=6+ to 8+ and 10+ states 
in lo^Sn. 

4.6 Tin isotopes for 100 < < 132 

In Fig. 1171 we present the model space and pertinent single-particle energies for a shell-model space for neutron valence 
holes using ^'^^Sn as closed shell core. In Table|5]we display selected states. As can be seen, the well-known near constant 
0+ — 2+ spacing is well reproduced, all the way down to ^^®Sn. Also the additional calculated states are in very good 
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Q-space: TV > 82 




Figure 17: Shell-model space for tin isotopes for 100 < A < 132. The valence neutrons are hole states. 



agreement with experiment. However more detailed analysis of the results close to ^^^Sn indicates that our effective 
two-particle interaction has difficulties in reproducing the shell closure which is believed to occur in this region. The 
increase of the the 0"*" — 2"*" splitting is not as sharp as found experimentally, even if the phenomenon is rather weak in 
the case of Sn. We have observed a similar feature around "^^Ca JHI which is generally agreed to be a good closed shell 
nucleus. There the deviation between theory and experiment is severe. Preliminary analysis indicates that our effective 
interaction may be slightly too attractive when the two particles occupy different single-particle orbits. This may be 
related to the radial wave functions which in our calculation are chosen to be harmonic oscillator functions. Further 
analysis of these spectra can be found in Refs. "99' . 

It is also instructive to study the degree of pairing in these systems. We calculate the squared overlaps between the 
constructed generalized seniority states in Eqs. H107|l and (jlOSfl and our shell model states through 



(f-2) 



|(^Sn(SM);0+|(5t)t|0)|2, 
K^Sn(SM);J,|i?t^,(5t)f-i|o) 



(120) 



The vacuum state |0) is the -'^'^^Sn-core and n is the number of valence particles. These quantities tell to what extent the 
shell model states satisfy the pairing picture, or in other words, how well is generalized seniority conserved as a quantum 
number. 

The squared overlaps are tabulated in Tabled and vary generally from 0.95 to 0.75. As the number of valence particles 
increases the squared overlaps are gradually decreasing. The overlaps involving the 4+ states show a fragmentation. In 
^^*Sn, the 4^ (SM) state is mainly a seniority v = 2 state. As approaching the middle of the shell, the next state, 4^, 
takes more and more over the structure of a seniority v = 2 state. The fragmentation of seniority over these two states 
can be understood from the fact that they are rather close in energy and therefore may have mixed structure. 



4.7 Tin isotopes above A = 132 

For tin isotopes above A = 132 calculations were performed for isotopes with 134 < A < 139 taking Z = 50, iV = 82 
as core. The active P-space for particles includes the I/7/2, lp3/2j 2pi/2, O/ig/2, I/5/2 and 0ii3/2 particle orbits. The 
single-particles energies e{fy^) = 0.00 MeV, £(^3/2) = 0.8537 MeV, e{py^) = 1.5609 MeV, eif^/^) = 1-6557 MeV 
e(/i~^2) = 2.0046 MeV and e(i^3/2) = 2.81 MeV are extracted from |10()j . The obtained results are presented in Table 01 
for isotopes with even values of A and in Table for isotopes with odd A. 

There are however no experimental data for the A > 134 tin isotopes available at the present time. We therefore 
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Table 2: Exitation spectra for the heavy Sn isotopes. 
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Table 3: Seniority i; = overlap | O+KS-t) 1 10) ^ and the seniority v = 2 overlaps \{'^Sn- J f\Dlj^j{S'^)'i-'^\Q)\'^ for 
the lowest-lying eigenstates of ^^^~^^°Sn. 
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cannot make conclusions on how good our effective interaction is for this mass region. A more detailed analysis of these 
results can be found in Ref. jl01| . 

4.8 Light antimony isotopes 

The nucleus ^°^Sb was produced in the reaction ^°Cr(^*Ni,lp2n) at a beam energy of 225 MeV with a 2.1 mg/cm^ thick 
target. The experiment was performed with the GAMMASPHERE Ge-detector array [TU^ at the ATLAS accelerator 
at Argonne National Laboratory. The experimental setup consisted of 78 Ge detectors, 95 Csl scintillators known as 
Microball Ji)3, for light charged particle detection, and the newly developed Neutron Shell. 

The calculation favors a 7^^=5/2+ assignment for the ground state in agreement with the suggestion from the proton 
decay data. In this state the valence proton is mainly in the orbit and the two neutron pairs are almost evenly 
distributed over the and g'j/2 neutron orbits. The situation is very similar in the 9/2+ and 13/2"*" states, while the 
z^dg^2g7/2 configuration exhausts the largest parts of the wave functions of the 15/2+ and 17/2+ states. The neutron part 
of the wave function of the 19/2+ state is almost identical to the 17/2+ state. However, since 17/2+ is the maximum spin 
for the ''^'^\/2^'^\/2^\ /2 configuration, the odd proton resides almost exclusively in the g7/2 orbit in the 19/2+ state. The 
19/2+ state is therefore the lowest lying of the observed states, in which the ng'j/2 orbit has a significant contribution to 
the wave function. This pattern repeats itself for the states with spins 21/2+, 23/2+, 25/2+ and 27/2+ in an alternating 
fashion with the state 21/2+ having a proton in the single-particle orbit TTd^/2- For proton degrees of freedom the Si/2, 
d3/2 and \iii/2 single-particle orbits give essentially negligible contributions to the wave functions and the energies of the 
excited states, as expected. For neutrons, although the single-particle distribution for a given state is also negligible, 
these orbits are important for a good describtion of the energy spectrum, as also demonstrated in large-scale shell-model 
calculations of tin isotopes [HHI ■ The effective interaction employed here is the same as that employed in Ref. [^Hl ) with 
an in general good reproduction of the data. 

We note here as well that the agreement with the experimentally proposed spin assignements is very good. The 
reason for such a good agreement is most likely that the wave functions of the states are dominated by neutronic degrees 
of freedom. The unbound proton is only a spectator while the well bound neutrons change orbits and alignment in 
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Table 4: Low-lying states for ^^^Sn, ^^^Sn and ^^*Sn. The energy eigenvalues are sorted according to the angular 
momentum assignment. 
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Figure 18: Proposed level scheme for ^°^Sb. Shell model calculations are shown on the right hand side. The widths of 
the arrows are proportional to the intensity of the transitions. 
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Table 5: Low-lying states for -'^^^Sn and ^'^'Sn. 
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transitions from high spin states. 

The recently described spectra of ^°^Sb show also low- lying 7/2+ and 11/2+ states, at approximately the same 
excitation energies as those shown in the theoretical calculation here. The 7-ray transitions feeding these two states in 
^"^Sb are as in ^°^Sb expected to be much weaker than the main 7-ray cascade, which explains why these two states were 
not identified in the experiment. They are interesting since their wave functions contain mainly contributions from the 
7rg7/2 orbit that is poorly known near ^""Sn, but a more sensitive experiment is needed for their identification. 

4.9 Neutron-rich nuclei in the IsOrf-lpO/ shells 

Studies of extremely neutron- rich nuclei have revealed a number of intriguing new phenomena. Two sets of these nuclei that 
have received particular attention are those with neutron number N in the vicinity of the IsOd and O/7/2 shell closures 
{N K, 20 and N k, 28). Experimental studies of neutron-rich Mg and Na isotopes indicate the onset of deformation, 
as well as the modification of the A'^ = 20 shell gap for ^^Mg and nearby nuclei |104| . Inspired by the rich set of 
phenomena occurring near the = 20 shell closure when N ^ Z, attention has been directed to nuclei near the = 28 
(sub)shell closure for a number of S and Ar isotopes |105l I106j where similar, but less dramatic, effects have been seen 
as well. In parallel with the experimental efforts, there have been several theoretical studies seeking to understand 



Table 6: The computed and measured values of B{E2) for the nuclei in this study using Cp = 1.5 and e,i = 0.5. 
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and, in some cases, predict properties of these unstable nuclei. Both mean-field |107[ I108j and shell-model calculations 
|1()5I llOfil 11091 [TTT)I lllll [il2l [ll3 have been proposed. The latter require a severe truncation to achieve tractable model 
spaces, since the successful description of these nuclei involves active nucleons in both the IsOd- and the IpO /-shells. The 
natural basis for the problem is therefore the full IsOd-lpOf space, which puts it out of reach of exact diagonalization on 
current hardware^. 

Shell-Model Monte Carlo (SMMC) methods [3 11141 I115| offer an alternative to direct diagonalization when the 
bases become very large. Though SMMC provides limited detailed spectroscopic information, it can predict, with good 
accuracy, overall nuclear properties such as masses, total strengths, strength distributions, and deformation, precisely 
those quantities probed by the recent experiments. 

There is limited experimental information about the highly unstable, neutron-rich nuclei under consideration. In many 
cases only the mass, excitation energy of the first excited state, the B{E2) to that state, and the /3-decay rate is known, 
and not even all of this information is available in some cases. From the measured B{E2), an estimate of the nuclear 
deformation parameter, (32, has been obtained via the usual relation 



with Ro = l.2A^/^ fm and B{E2) given in e^fm^. 

Much of the interest in the region stems from the unexpectedly large values of the deduced /32, results which suggest 
the onset of deformation and have led to speculations about the vanishing of the = 20 and N = 28 shell gaps. The 
lowering in energy of the 2^ state supports this interpretation. The most thoroughly studied case, and the one which most 
convincingly demonstrates these phenomena, is ■^^Mg with its extremely large B{E2) = 454 ± 78 e^fm'* and corresponding 
P2 = 0.513 |1()4| : however, a word of caution is necessary when deciding on the basis of this limited information that we 
are in the presence of well-deformed rotors: for ^^Mg, we would obtain P2 — 0.67, even more spectacular, and for ^^C, 
(32 = 0.8, well above the superdeformed bands. 

Most of the measured observables can be calculated within the SMMC framework. It is well known that in deformed 
nuclei the total B{E2) strength is almost saturated by the 0+, 2^ transition (typically 80% to 90% of the strength 
lies in this transition). Thus the total strength calculated by SMMC should only slightly overestimate the strength of 
the measured transition. In Table 1 the SMMC computed values of B{E2, total) are compared both to the experimental 
B{E2;0'^^ 2f) values and to the values found in various truncated shell-model calculations. Reasonable agreement 
with experimental data across the space is obtained when one chooses effective charges of = 1.5 and e„ = 0.5. All 
of the theoretical calculations require excitations to the lpO/-shell before reasonable values can be obtained. We note 
a general agreement among all calculations of the B{E2) for ''^Ar, although they are typically larger than experimental 
data would suggest. We also note a somewhat lower value of the B{E2) in this calculation as compared to experiment 
and other theoretical calculations in the case of ''^S. Table 2 gives selected occupation numbers for the nuclei considered. 
We first note a difficulty in extrapolating some of the occupations where the number of particles is nearly zero. This leads 
to a systematic error bar that we estimate at ±0.2 for all occupations shown, while the statistical error bar is quoted 
in the table. The extrapolations for occupation numbers were principally linear. Table 2 shows that ^^Mg remains as 
an almost pure sc?-shell nucleus, as expected. We also see that the protons in ^°Ne, ^^Mg, and ''^Si are almost entirely 
confined to the sd shell. This latter is a pleasing result in at least two regards. First, it shows that the interaction does 
not mix the two shells to an unrealistically large extent. Second, if spurious CoM contamination were a severe problem, 
we would expect to see a larger proton O/7/2 population for these nuclei due to the Oc?5/2-0/7/2 "transition" mediated by 
the center-of-mass creation operator. The fact that there is little proton f'j/2 occupation for these nuclei confirms that 
the CoM contamination is under reasonable control. See Ref. for further details. An interesting feature of Table 2 
lies in the neutron occupations of the iV = 20 nuclei (^''Ne and '"'^Mg) and the TV = 28 nuclei (''^Si, "45, and ^^^Ar). The 
neutron occupations of the two = 20 nuclei are quite similar, confirming the finding of Fukunishi et al. [lllj and Poves 
and Retamosa |110| that the = 20 shell gap is modified. In fact, the neutron O/7/2 orbital contains approximately 
two particles before the TV = 20 closure, thus behaving like an intruder single-particle state. Furthermore, we see that 
2p-2h excitations dominate although higher excitations also play some role. We also see that the neutrons occupying the 
IpO /-shell in iV = 20 systems are principally confined to the O/7/2 sub-shell. The conclusions that follow from looking at 
nuclei with N > 20, particularly those with N = 28, are that the = 20 shell is nearly completely closed at this point, 
and that the N = 28 closure shell is reasonably robust, although approximately one neutron occupies the upper part of 
the IpOf shell. Coupling of the protons with the low-lying neutron excitations probably accounts for the relatively large 
B{E2), without the need of invoking rotational behavior. In Table 3 we show the SMMC total Gamow- Teller (GT~) 
strength. We compare our results to those of previous truncated calculations, where available. In all cases, our results 
are slightly smaller than, but in good accord with, other calculations. Since we do not calculate the strength function, 
we do not compute /3-decay lifetimes. 

®For a treatment of the CoM problem, see Ref. 1691 . 
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Table 7: The calculated SMMC neutron and proton occupation numbers for the sd shell, the O/7/2 sub-shell, and the 
remaining orbitals of the pf shell. The statistical errors are given for linear extrapolations. A systematic error of ±0.2 
should also be included. The first row represents neutron results, while the second row represents protons. 





N, Z 


IsOd 


O./7/2 


lpO/5/2 


— OOti 1 

22]V[g 


10,12 


3.93 ± 0.02 


0.1 ± 0.02 


-0.05 ±0.01 






2.04 ± 0.02 


0.00 ± 0.01 


-0.05 ±0.01 




20,10 


9.95 ± 0.03 


2.32 ± 0.03 


-0.26 ±0.02 






2.03 ± 0.02 


-0.01 ± 0.01 


-0.02 ±0.01 


32]VIg 


20,12 


9.84 ± 0.03 


2.37 ± 0.03 


-0.21 ± 0.02 






3.99 ± 0.03 


0.05 ± 0.02 


-0.05 ±0.01 


36 Ar 


18,18 


9.07 ± 0.03 


1.08 ± 0.02 


-0.15 ±0.02 






9.07 ± 0.03 


1.08 ± 0.02 


-0.15 ±0.02 


4"S 


24,16 


11.00 ± 0.03 


5.00 ± 0.03 


-0.01 ± 0.02 






1.0 1 ± 0.04 


0.54 ± O.Oz 


A 1 _L A AO 

— O.lz ± O.Oz 


42 Si 


28,14 


11.77 ±0.02 


7.34 ±0.02 


0.90 ±0.03 






5.79 ±0.03 


0.25 ±0.02 


-0.07 ±0.01 


42 s 


26,16 


11.41 ±0.02 


6.33 ±0.02 


0.25 ±0.03 






7.49 ±0.03 


0.58 ±0.02 


-0.09 ±0.02 


44s 


28,16 


11.74 ±0.02 


7.18 ±0.02 


1.06 ±0.03 






7.54 ±0.03 


0.56 ±0.02 


-0.12 ±0.02 


44Ti 


22,22 


10.42 ±0.03 


3.58 ±0.02 


0.00 ±0.02 






10.42 ±0.03 


3.58 ±0.02 


0.00 ±0.02 


46 Ar 


28,18 


11.64 ±0.02 


7.13 ±0.02 


1.23 ±0.03 






8.74 ±0.03 


1.34 ±0.02 


-0.08 ±0.02 



Table 8: The calculated total Gamow- Teller strength, GT , from this study. The results of other studies, when available, 
are presented for comparison. 



Nucleus 


SMMC 


Other 


2^Mg 


0.578 ±0.06 




30Ne 


29.41 ±0.25 




32 Mg 


24.00 ±0.34 




36 Ar 


2.13 ±0.61 




40s 


22.19 ±0.44 


22.87 112 


42 s 


28.13 ±0.42 


28.89 ^112 


42Si 


40.61 ±0.34 




44s 


34.59 ±0.39 


34.93 .112 


44Ti 


4.64 ±0.66 




46 Ar 


29.07 ±0.44 


28.84 112 
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5 Non-perturbative resummations: parquet diagrams 



5.1 Inclusion of hole-hole contributions and single-particle propagators 

With the G-matrix defined according to the double-partitioned scheme we can easily solve Eq. I)29f) through matrix 
inversion. The number of hole-hole and particle-particle configurations is then rather small, typically smaller than ~ 100, 
and a matrix inversion is then rather trivial. Before we discuss the solution of Eq. H29I) , it is always instructive to consider 
the contributions to second order in perturbation theory, i.e., diagrams (a) and (b) of Fig.[Sl The external legs can be 
particle states or hole states. Diagram (a) reads 

(122) 

pq PI 

and 

a/3 

We note here the minus sign in the energy denominator, since in the latter expression we are using the hole-hole term of 
the propagator of Eq. (|30|l . If we use a double-partitioned G-matrix for say ^^O and are interested in an effective valence 
space interaction for the IsOd-shell^, then typically the single-particle orbits of the intermediate states will be represented 
by states in the IpOf major shell. Hole states are then defined by single-particle states in the Os and Op shells. Clearly, 
the number of two-body intermediate states is rather limited. To third order we have diagrams like (c) and (d) of Fig.[S] 
Diagram (c) is just the third-order equivalent of Eq. 1)1221) and reads 

- A 2^ ^i'ipqJ ^pqrwJ „ _ ^ _ ^ ^rwUJ^ U^^i 

pqrw ^ ^ 

while diagram (d) contains both a two-particle and a two-hole intermediate state and reads 

Thus, solving Eq. (|29|l will then yield contributions to the effective interaction such as the above expressions. 

Here we have also tacitly assumed that the energy denominators do not diverge, i.e., we have chosen an energy s so 
that we avoid the poles. This has always been the standard approach in calculations of shell-model effective interactions. 
To give an example, consider now diagram (b) and suppose that we are using harmonic oscillator wave functions. Let us 
also assume that the two hole states are from the Op-shell and that the valence particles are in the IsOd-shell. If we rescale 
the energies of the valence space to zero, then the two-hole state would yield —28 MeV with an oscillator parameter 
h = 1.72 fm. If s = —28, the denominator diverges. In this case it is rather easy to obtain the imaginary part, and even 
if we were to chose s different from —28 MeV, the imaginary part will influence the real part of the effective interaction 
through dispersion relations, see e.g., Refs. [491 15UL ri20) . It is therefore at best just a first approximation to neglect the 
imaginary term. Moreover, if we solve Dyson's equation for the self-energy, the single-particle energies may contain an 
imaginary part. Technically it is however not difficult to deal with imaginary contributions, one needs to invert a complex 
matrix rather than a real one. These technicalities will however be described elsewhere |121| . 

Using the double-partitioned G-matrix, we can then rewrite Eq. (|29|l as 

fs^-G^^^i +iVG[^^i c[i2]r[i2] Q26) 

56 

where G'^^I jg j^gt the double-partitioned G-matrix discussed above. It is also energy dependent, in contrast to V . In case 
we were to employ this equation for effective interactions in the IsOd-shell, the intermediate two-particle states would then 
come from just e.g., the IpO /-shell. This equation, which now is solved within a much smaller space than the original 
one spanned by the total Qpp, allows clearly for computationally amenable solutions. It corresponds to the so-called 
model-space approach to the solution of the Feynman-Galitskii equations. Thus, a possible approach would consist of the 
following steps 

1. Solve the G-matrix equation from Eq. (|50(l using the double-partitioning scheme. 

2. The next step is then to solve Eq. 129(1 and Dyson's equation for the self-energy. 

3. This scheme is iterated till self-consistency is achieved, see the discussion below. 



''This means that the labels 1234 will refer to particle states in the IsOd-shell. 
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We will however not employ this model-space scheme in our actual calculations. There are several reasons for not 
doing so. 

Let us first assume that we omit the [f3] and [14] channels in our iterative scheme for Eq. 1)126(1 . The next iteration 
of Eq. 1(126(1 would then look like 

^[12] _ -p[12] p[12]^[i2]p[12] „„x 

(1) ^ (0) + ^ (0) ^ (1) ' y^"^' ) 

[l2l *TTT?Tli [l2l 

where the vertex function F^q^' is the solution of Eq. 1(126(1 . However, we cannot define the "bare" vertex F^p^' to be 

the solution of Eq. I(126() simply because then we would be double- counting contributions. Thus, F^g^' has to equal the 
G-matrix. The only change in Eq. ((127() arises from the solution of Dyson's equation and thereby new single-particle 
energies. 

Let us then for the sake of simplicity assume that the single-particle energies are just the Hartree-Fock solutions. 
The problem we are aiming at arises at the Hartree-Fock level. In order to obtain Hartree-Fock solutions which are 
independent of the chosen harmonic oscillator parameter b, we typically need to include single-particle orbits from quite 
many major shells. Typical constraints we have found when we do so-called Brueckner-Hartree-Fock (BHF) calculations 
for finite nuclei is that we need at least 2n + I < 10 in order to obtain a result which is independent of the chosen 
value. The way we solve the BHF equations is to expand the new single-particle wave functions tpx, with A representing 
the quantum numbers nlj , in terms of harmonic oscillator wave functions, i.e., 

2n+l 

IV'a) = E Ci'^ l'^") (128) 

where (f)a are the harmonic oscillator wave functions with quantum numbers a — nlj and C are the coefficients to be 
varied in the Hartree-Fock calculations. The single-particle energies at the Hartree-Fock level are just 

-f^(a/i|G(£„-l-e,,)|a/i), (129) 

h 

where the single-particle states are just those of the harmonic oscillator. The G-matrix used in the first iteration in the 
BHF calculation is the one given by the solution of Eq. ((50(1 . The coefficients Ca can then be obtained by diagonalizing 
a matrix of dimension N x N, where N is the number of single-particle orbits with the same Ij values. 



5.2 Screening corrections and vertex renormalization, the equations for the [13] and [14] channels 

We start as in the previous section with the definition of the interaction vertices in the [13] and [14] channels and the 
corresponding integral equations. Thereafter, we discuss various approximations to these equations such as the summation 
of TDA and RPA diagrams. Eventually, the aim is to merge the discussion in this section and the preceeding one into 
equations for a self-consistent scheme which combines all three channels, namely the so-called set of parquet equations to 
be discussed in sectional 

The equations for the renormalized vertex in the [13] and [14] channels have the same form as Eq. ((28(1 . namely 

r[13] ^^[13]^^[13](^^)p[13]^ (^30) 

and 

fM ^V^[i4] + ^[i4](^^^)r[i4]. (131) 

The matrix elements which enter are however defined differently and the irreducible diagrams of T^I^'^l and T^I^^l can 
obviously not be the same. With irreducible in the [13] channel we will mean a diagram, which by cutting an internal 
particle-hole pair, cannot be separated into a piece containing the external legs 1,3 and another piece containing 2,4 as 
external legs. The definition for the irreducible vertex in the [14] channel is similar and we illustrate these differences in 
Fig. 1191 Diagram (a) is just the lowest-order interaction in the [13] channel and is therefore irreducible. Diagram (b) is 
an irreducible diagram in the [13] channel, whereas it is reducible in the [14] channel. Diagram (c) is in turn irreducible 
in the [14] channel and reducible in the [13] channel. Diagram (d) is an example of a diagram which is irreducible in both 
channels. This diagram stems from the [12] channel. 

The energy variables in these channels are, following Fig. O and Eqs. 1(17(1 and ((18(1 . 

t = £3 - £1 = £2 - £4, (132) 

for the [13] channel and 

u = £1 - £4 = £3 - £2, (133) 

^Throughout this work our unperturbed single-particle basis par excellence will always be that of the harmonic oscillator. 
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(a) (b) (c) (d) 

Figure 19: Examples of irreducible and reducible diagrams in the [13] and [14] channels. See text for further details. 



for the [14] channel. Defining the unperturbed particle-hole propagators in the energy representation as 
and 



' hp 



t - Ep + Eh + tr] t + Ep - Eh 



IT] 



;[14] 



u- Ep + Eh + ir] u + Ep - Eh - trj 
we arrive at the following equations for the interaction vertex in these two channels 

-[13] _ T/[13] ^Y^T/[13] ^[13]p[13] 



r 



1234.J\'') ^ ^1234,7 ' ^ I2phj^ 

ph 



p/i34 



and 



.[14] 
1234J 



{u) = V, 



[14] 
1234,7 



/ . *12phJ- 
ph 



pMAJ 



in). 



(134) 



(135) 



(136) 



(137) 



These equations, together with Eq. (|29|l . can then form the basis for the first iteration in a self-consistent scheme for 
renormalization corrections of the parquet type. The origin of the minus sign in Eq. H137() follows from the diagram rules 
and will be examplified below. A graphical view of these equations is given in Fig. [201 The reader should also keep 
in mind the two contributions to the particle propagators of Eqs. H134() and H135() . 



5.3 Screened ph and 2p2h interactions 

Here we study the screening of the particle- hole and the 2p2h interactions given in Fig. ^ indicated by Vph and V2p2h, 
respectively. Before we list the final expression, it is however instructive to consider the corrections to second order in the 
interaction V to the ph and 2p2h vertices. In Fig. |^ we display the second-order corrections to the ph diagrams of Fig. 
n Diagram (a) is the core-polarization correction term to the particle-hole interaction, and corresponds to a contribution 
from the [14] channel, as indicated by the coupling order. The term labeled (b) corresponds to the exchange term of (a) 
and is coupled in the [13]-order, see also the discussion in connection with Eqs. 1^511 and (|^ . The other corrections, like 
(c) and (d) include particle-particle and hole-hole intermediate states, respectively. They are irreducible in both the [13] 
channel and the [14] channel, and can therefore enter the irreducible vertices of these two channels in later iterations. 
They are however not generated by various iterations of Eqs. (|136|l and H137|l . In fact, if we replace F by G in Eqs. H136(l 
and (|137|) , diagram (c) is already accounted for by the G- matrix. It may however be included if the double-partitioned 
G-matrix of the previous section is used. Let us now look at the analytical expressions in an angular momentum coupled 
basis for diagrams (a) and (b) of Fig. [2] Here we just include the first term of the propagators of Eqs. (|134|l and (|135|l . 
The second terms will give rise to the 2p2h contributions discussed below. In the following discussion we will also assume 
that the interaction V does not depend on the energy, although it is rather easy to generalize to an energy dependent 
interaction. Diagram (a) reads 

(«) = - E(-)"^^-'(-)^^-d./ , + / _ y^Zj^ (138) 
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Figure 20: (a) shows the structure of the integral equation for the interaction vertex in the [13] channel, (b) represents 
the integral channel for the [14] channel. The coupling order is displayed as well. 
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(d) 



Figure 21: Second-order perturbation theory corrections to the ph interaction vertex. 
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The factor (— )^-'^ stems from the opening up and recoupUng of an internal particle-hole pair and the phase 
is needed in order to rewrite the matrix elements in the coupling order of Eq. 1)17(1 . The general structure of Eq. H138|l is 
just of the form -K'^'Qi^Ve'^'K'"' , with eli^l = £^ + - Sp ^ Er ^ u + . 



Er and we have defined 



(139) 



for the on-shell energy case. This is the equivalent of the energy variable of Eq. 127|l in the [12] channel. Diagram (b) is 
in turn given by 

(140) 



r7 



and we note that the contributions are clearly different. The minus sign in Eq. p38|) stems from the standard diagram 
rules [SSj. In our use of the diagram rules below, we will omit the use of the rule for the number of external valence hole 
lines. In our case then, as can also be deduced from inspection of Fig. |^ diagram (a) has zero closed loops and three 
hole lines, giving thereby rise to a minus sign. Diagram (b) has an additional closed loop and thereby yielding the plus 
sign. The energy denominator is in this case 

e[i3l =t + £^-er, (141) 

with 

t = Eg - Ep ^ Ep- Ea- (142) 

We notice, using the relations discussed in Eqs. (|23|l and H24|l . that diagram (a) is simply the exchange diagram of (b). 
We need however to include both diagrams in order to obtain an antisymmetric equation for the particle-hole channels 
which exhibits the same properties as the [12] channel shown in Eq. H26|l . This is actually crucial in solving the parquet 
equations. We wish namely that every iteration, with a given approximation to the vertex function V , preserves the 
antisymmetry property. This point cannot be emphasized enough. Let us now see what happens to third order in the 
interaction. Third order corrections to the ph vertices (a) and (b) involving only ph intermediate states are shown in (a) 





(a) 



(b) 



Figure 22: Corrections beyond second order in the interaction V to the ph interaction vertex, (a) is in the [14] channel 
and (b) is in [13] channel. 



and (b) of Fig. 1221 respectively. The analytical expression for the third-order contribution (a) is given by 



rs7(5 



U + Ey ~ Er 



[14] 
S'yJ 



U + ES - Es 



[14] 



sPqSJ ^ 



(143) 



with / — (^—yr+jB+j-y+js 2J ^_yj.y+2j5 ^ rpj^jg equation has the general structure 

, .[14] ^Vph_ [14] Vph_ [14] 
ph g[14] "^Ph g[14] '^ph • 

A similar expression applies to diagram (b), whose expression is 



rs^S 



[13] 
* psq&J- 



(144) 
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It has the general structure 

^^[13] '^ph ,^[13] '■^ph ,^[13] 
>h g[13] '^ph g[13] '^ph ■ 

But these expressions have the same sign! Diagram (a) counts now 4 hole lines, and (b) counts also 4 hole lines and 2 
closed loops. However, there are three interaction terms V , and taking the exchange term of each of these in diagram (a) 
leads to the desired results, namely (o) = — (&), as it should. Thus, to third order we keep the antisymmetry property of 
r in the [13] and [14] channels. It is easy to see that in the [14] channel we will always have an alternating sign in front 
of each contribution, since every new order in perturbation theory brings a new hole line and no closed loop, and thus a 
new minus sign. In the [13] channel we have always one new hole line and one new closed loop for every new vertex. If 
we consider only the screening of the ph vertex, we can then set up a pcrturbativc expansion in terms of the ph vertex 
for the vertex functions F^^'^l and pl^^l. For notational economy we will skip the Pauli operators Qp^ ^p discussions 
below. It will always be understood that the intermediate states are two-body particle-hole states, |ph) or |hp). Consider 

e.g., r["i 



^ - ^'ph g[i4] '^ph + Vh ^[14] g[i4] '^ph + • • ■ ' U^OJ 

which can be summed up to yield 



r[14] _ y[14] _ y[14] 1 T.[14] _ y[14] _ ^[14] ^p[14] Q4g^, 

which is the standard TDA expression for the ph term. The corresponding expression in the [13] channel results in 

r[^^i = C^ + C^r[-i. (147) 

The signs agree with the expressions of Blaizot and Ripka see chapter 15 and Eq. (15.50). The summations in both 
channels ensures that the final vertex is antisymmetric and the combination of the latter two equations results in the 
familiar TDA equations. We next look at the 2p2h matrix element and show the corresponding corrections to second 
order in perturbation theory in Fig. 1231 If we omit diagrams (e) and (f) which contain 2p and 2h intermediate states 



(a) 



P q 



(b) 







(d) (e) (f) 

Figure 23: Corrections to second order in V of the 2p2h vertex. 



generated by the solutions in [12] channel, we have for diagram (a) 



[14] 



1i I ^ T 



with the general structure 



[14] J_T/[14] 



2p2h ^[14] *^2p2h 



(148) 



(149) 



44 



Note well the minus sign in front of u. The contribution from the propagator can in this case be retraced to the second 
term of the propagator of Eq. H135|l . Diagram (a) follows the coupling order of the [14] channel. It is also easy to see that 
diagram (b) is given by 

ib) = -E(-)""^^--'(-)'^-^'X/ _„^! (150) 



-U + 

r7 ' 



and has the structure 



V g[14] %2h- UO^j 



Similarly, if we now move to the [13] channel we have the following expressions 



r7 

and 



r7 

with the general structure 
and 



id) = E(-)"^^'^"'(-)'^'^Cz7-f-r^S.' (153) 



yim _J_y[13] 

'^2p2hg[13] ''2p2h' (^O^) 
V g[13] %2h' UOOj 

respectively. Diagram (c) is just the exchange of (a) and includes two 2p2h vertices, while diagram (d) is the exchange of 
diagram (b) and includes a ph vertex multiplied with a 2p2h vertex. We note again that the antisymmetry is ensured at 
a given order in the interaction only if we include the corrections at the same level in both channels. 

One can then easily sum up higher-order corrections to the 2p2h diagrams as well in both channels. The inclusion of 
the backward going particle- hole pair in the propagators of Eqs. (|134() and H135|l ensures thus that we will also sum to 
infinite order 2p2h corrections. This leads ultimately to the familiar RPA equations, see e.g., Refs. |48|. 

A closer inspection of Eqs. (|15U|I and 1)1531) shows that if we only include ph vertices, we could resum these corrections 
to infinite order for the 2p2h vertex by observing that the structure of such diagrams would be of the form (e.g., in the 
[14] channel ) 

p[14] _ T.[14] _ T^[14] J_^[14] ^[14] J_^[14] J_^[14] 

^ 2D2h — ''2D2h "oh ,f1 41 '^2D2h ^ "^oh ,fl4l''^Dh fl 41 '^2D2h ' ' ' ' ' \^^^) 



which can be summed up to yield 

^[14] _ [14] [U] 1 [14] , . 

2p2h - ''2p2h g[i4] _ ^[14] *^2p2h' l^O/j 

and similarly for the [13] channel, but with a plus sign. The modification discussed in Eqs. H146|l and H157|l serve to 
modify the propagation of a particle- hole pair and have normally been termed for propagator renormalizations, as can 
easily be seen from Eqs. (|146|) and H157|) where the propagation of a free particle-hole pair is modified by the presence of 
the interaction V in the energy denominator. Other important processes which can affect e.g., various polarization terms 
are those represented by so-called vertex renormalizations, a term originally introduced by Kirson and Zamick |122| . 
These authors studied the renormalizations of the 2plh and 2hlp vertices as well, see also Kirson 07] and Ellis and Osnes 
0H| for further discussions. We will therefore end the discussion in this section by looking at such renormalizations. 



5.4 Further renormalizations 

In the previous subsection we dealt mainly with what has conventionally been labelled for propagator renormalizations. 
We will therefore extend the standard TDA and RPA scheme by looking at further ways of renormalizing interaction 
vertices. The approach discussed here follows Kirson 47 . Extensions were made later by Ellis and Goodin |123| and 
Ellis, Mavromatis and Miither We will limit the discussion here to the scheme of Kirson. We start therefore with 
the contributions to second order to the 2plh vertex^. These contributions are shown in Fig.[^ Diagram (a) consists of 
a 2plh vertex multiplied with a ph vertex whereas (b) stems from the multiplication of a 2p2h vertex with a 2plh vertex. 
They both contain a particle-hole pair as an intermediate state and follow the coupling order of the [14] channel. The 
exchange diagram of (b) is given by (c), while that of (a) is diagram (d). Diagrams (e) and (f) represent contributions 
from the [12] channel and are irreducible in both particle-hole channels. Before we sketch the general structure of the 

®The discussion here applies to the other interaction vertices discussed in Fig. [Tl 
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Figure 24: The corrections to second order in V of the 2plh vertex. 



renormalization procedure of Kirson, it is instructive to consider again the equations to second order in perturbation 
theory, as the general expressions can be deduced from inspection of these diagrams. Diagram (a) is 

(a) - ^(-)^--+..-./(_)2,-.^[i41 1 ^[14]^^ (,53) 

and the plus sign stems from the diagram rules |58| . i.e., we have two hole lines and no closed loop. The propagator is 
that arising from the first term in Eq. H135|) . The general structure is 

Diagram (b) reads 

ib) = E(-)"^^-^'(-)'^-C _„^! (160) 



with the following structure 



J_y[i4] ^gl^ 

^2plhg[14] ''2p2h- (^'^^) 

In this case the propagator stems from the second term in Eq. (|135|) . Similar equations arise for e.g., the 2hlp vertices of 

Fig.m 

Before we write down the self-consistent equations of Kirson 47 , let us assume that we can approximate the 2plh 
vertex in the [14] channel by the first order term and diagrams (a) and (b). This renormalized vertex, which we here 
label V2pih, is then given by 

V2plh ~ V2plh + Vh ^[n]-'^2plh + %lhgnf '^2p21i- Uo^j 

If we now allow for the screening to infinite order of the ph vertex given by Eq. Ij 146(1 and replace the 2p2h vertex in the 
above equation with Eq. ((157|l we obtain the following renormalization of the 2plh vertex 



t7 _ ta[14] I T^[14] 4 [14] [14] 1 [14] 

V2plh - V2plh + V 414] _ r^[14] %lh + ^2plh _ [14] ^^2p2h- i^^dj 



A similar equation applies to the 2hlp vertex of Fig. ^ and for the [13] channel. Eqs. (|146|l . H157() and (|163|l form then 
the starting point for the approach of Kirson |47| . Examples of diagrams which can be obtained through the iterative 
solution of Eqs. (|146|l . (|157|l and (|163|l are given in Fig. [23 

The question now is however how to relate Eqs. H146|) . H157|l and (|163|) with those from Eqs. p36|) and (|137|l . This is 
rather trivial if we recall that the labels 1234 can, as was also discussed in section |2 represent whatever single- particle 
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(a) 



(b) 



Figure 25: Examples of diagrams which can arise from Kirson's self-consistent equations. 



states, either holes or particles. Thus, ^1234 can represent a 2plh, 2hlp, 2p2h, 2p, 2h or a ph vertex. This means that, due 
to the choice of propagators in Eqs. H134|) and (|135|l equations like Eq. Hlt)3|) are already inherent in Eqs. (|136|l and (|137|l . 
If we e.g., approximate Eq. H137|l to second order in the interaction V and let the single-particle labels 1234 represent a 
2plh interaction vertex, we immediately reobtain Eq. (|162|l . If we let 1234 represent a 2p2h vertex, we find to second 
order diagrams (a)-(d) of Fig. 1231 

Till now we have however refrained from discussing the contributions from the [12] channel, examples were only shown 
in (c) and (d) of Fig. |^ (e) and (f ) of Figs. |2S1 and These diagrams cannot be generated by simply iterating the 
equations for the [13] and [14] channels, but they could enter as contributions to the irreducible vertex in the [13] and 
[14] channels from the first iteration in the [12]-channel. We see thus the emerging contour of an iterative scheme. The 
crucial point is however how to perform the next iteration of say Eqs. I|146|) . H157|l and (|163|l or Eqs. 1)13(1 1) . H137|) and 
(I29|l from the [12] channel The question is how do we include the results from the first iteration into the next one, i.e., 
how to modify the bare vertices l/I^^l and T^I^'^l in e.g., Eqs. (|136|) and H137|) in order to obtain an effective interaction for 
the shell model. We have also not addressed how to deal with the solution of Dyson's equation for the one-body Green's 
function. We mention also that Ellis and Goodin 123 included pp correlations, i.e., terms from the [12] channel such 
as diagram (e) of Fig. |^ as well when they considered the screening of the 2plh and 2hlp vertices. Furthermore, as 
already mentioned in the introduction the authors of Refs. extended the pp RPA to include the particle-hole (ph) 

RPA, though screening of the 2plh and 2hlp vertices was not included. In Ref. [^3 however, a study with self-consistent 
single-particle energies was also performed. These works represent thus a first serious step towards the solution of the 
parquet equations, i.e., a many-body scheme which solves self-consistently the equations in the [12], [13] and [14] channels, 
with the addition of the self-consistent evaluation of the self-energy. It ought also to be mentioned that one of the really 
first applications for nuclear systems was performed in a series of papers by Dickhoff and Miither and co-workers |124| 
for nuclear matter. These authors actually performed the first iteration of the three channels. 

5.5 Parquet diagrams 

The equations we discussed in the two previous sections can be generalized in matrix form as 

r = r[*^i +r['^'i^[*^'ir, (i64) 

where obviously [ij] represents a given channel, Q^^^^ is the particle-particle, or hole-hole or particle-hole propagator. The 
propagator is a product of two single-particle propagators g which we do not specify any further here. They are defined 
by the solution of Dyson's equation in Eq. H31|l . The irreducible vertices must appear in the solution of the self-energy, 
and conversely, the self-energy must appear in all single-particle propagators within the expressions for the three channels 
[12], [13] and [14]. 

For all of our practical purposes, the irreducible vertex used in all channels is the G-matrix defined by the large model 
space of Fig. El 
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Let us now define the contribution from the [12], [13] and [14] channels following Ref. ^Uj. Eq. H29|l is then rewritten 

as 

L = r[i2l _ G, (165) 
where obviously L stands for ladder. The ladder term can then be rewritten as 

L = Gg^^^^G + Gg^^^^L. (166) 

In a similar way we can define the diagrams beyond first order in the particle-hole channel as 

^[13] ^ p[13] _ (lg7) 

and 

r[u] ^ r[i4] _ (^gg) 

where G now is coupled in either the [13] or [14] way. Rewriting i?^'-'! , where R refers to ring diagrams, we obtain 

rM ^ G^fe] G + G^l^^l , (169) 

where ij stands for either [13] or [14]. The equation for the vertex function F in Eq. (|164|l becomes then 

V = G + L + + i?!!"*! . (170) 

Note here that the vertex T can he represented in the coupling order of any of the above channels. Our convention is that 
of the [12] channel. 

If Eqs. p66|l and H169|) define the first iteration, it should be fairly obvious to see that the next iteration would be 

L={g + + i?[14]) gll2] ^ ^[13] ^ ^[14]^ ^[13] ^ ^[14]^ gll2]^^ (^7^) 

where the contributions R^^^^ + are recoupled according to the coupling order of the 12-channel. These contributions 
are irreducible in the [12] channel. Similarly, for the rings we have 

^(g + L + ^[13] (g + L + j + (g + L + g[i3]^[i4]^ (172) 

and 

i?M = (g + L + g[i4] (g + L + + (g + L + ^[i4]i?[i4]. (173) 

Our scheme for calculating T will be an iterative one based on Eqs. p7U|) - (|173|) and the solution of Dyson's equation 
for the single-particle propagator. This set of equations will then yield the two-body parquet diagrams. Relating the 
above equations to the discussions above it is rather easy to see that the G-matrix, TDA, RPA and Kirson's screening 
scheme are contained in Eqs. I|17U|) - (|173|I . Applications of this method will be presented elsewhere |121| . However, we 
present here a numerically viable approach to the parquet equations. Here we will limit oursclf to just sketch the structure 
of the solution, more technical details can be found in Ref. 1121). 

The iterative scheme starts with the solution of Ec^. (|29l) . We define then Q^^^l = Qpp + Qhh which results in 

The subscript (0) means that this is just the first iteration. The single-particle energies are the unperturbed harmonic 
oscillator energies. 

The first step in our calculations is to evaluate G. It is calculated following the discussion of sectional employing 
a large model space. The exclusion operator for this G-matrix was defined in Fig. El and is typically defined for 6—10 
major oscillator shells, as discussed in connection with Fig. There, 8 major oscillator shells were needed in order to 
obtain a converged result. Similar results are obtained in the coupled cluster calculations discussed in the next section. 

The second step is to solve Eq. ()174|l . which is now a complex equation. The hole states and particle states are then 
defined according to Fig. ^| These single-particle states define then the projection operator Q^^^l — Qpp + Qhh- This 
equation is then solved by matrix diagonalization, yielding a Green's function and the pertinent new vertex. 

These two steps lead then to the first iteration of the ladders, i.e., 

L(o) - T^l^} - G. (175) 

It contains both hole-hole and particle-particle intermediate states and is a complex matrix. The external single-particle 
legs can be particles or holes. Only unperturbed single-particle energies enter the definition of the two-body propagators. 
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The third step is to calculate the first iteration for the rings, namely 

R^^; = (G + L(o)) ^[1^1 (G + L(o)) + (G + L(o)) G^^'^R^i,';, (176) 

and 

i?M ^ (G + ^11^1 (G + L(„)) + (G + L(„)) ^I^^li^^ . (177) 

The equations for L and R are all defined within a truncated Hilbert space. They can therefore be recast into matrix 
equations of finite dimensionality. Recall also that we need to recouple the contribution from the [12] into the relevant 
ones for the [13] and [14] channels. This is done employing Eqs. (|17|l and (|18f) . With these contributions, we can now 
obtain the vertex function T after the first interaction 

r(o)=G + i(o)+i?g'+i?[o^ (178) 

The fourth step is to compute the self-energy and thereby obtain new single-particle energies. In so doing, care has to 
be exercised in order to avoid double-counting problems. A thourough discussion of this topic can be found in Ref. |4(Jj . 
More details will also be presented in Ref. jl21j . The new single-particle wave functions are obtained by diagonalizing a 
matrix of dimension Ua x Ua, Ua the quantum number n of the single-particle state a. 

The fifth step is to repeat steps 1-4 with the new single-particle energies till a predetermined self-consistency is 
obtained. But now the rings have to be included in all equations, i.e., we solve Eqs. H17UI) - (|173|I . 

The final vertex F can then be used to define a new effective interaction to be applied in shell model studies, where 
many more diagrams are considered than in present state of the art calculations, see e.g.. Fig. 8 of Ref. 40 for a list of 
diagrams to sixth order entering the definition of the irreducible vertex F. 

The same technique employing a large model-space is used in our coupled cluster calculations discussed in the next 
section. 

6 Non-perturbative resummations: Coupled cluster theory 

As we have seen from the above discussions, nuclear many-body theory often begins with a G-matrix interaction which 
is derived from an underlying bare nucleon-nuclcon interaction. In this section we limit the discussion to the no-core 
G-matrix so that all particles are active within our chosen model space. Using a given basis expansion of the many-body 
wave function we could then solve the nuclear problem by diagonalization as has been pursued by the No-Core shell 
model collaboration In fact, the current and most advanced no-core oscillator expansion techniques have 

approached ^^C, with nearly converged solutions ^125^ . 

It should be evident, however, that diagonalization procedures scale almost combinatorally with the number of particles 
in a given number of single-particle orbitals. Because of this scaling, diagonalization simply becomes untenable at some 
point. The efforts to expand diagonalization into p-shell nuclei with all nucleons active, an effort that spans over ten 
years, illustrates the problem. The computational complexity of the nucleus grows dramatically as the size of the nucleus 
increases. As a simple example consider oscillator single-particle states, and single-particle spaces consisting of 4 and 

7 major oscillator shells, and compare the number of uncoupled many-body basis states there are for 4,8,12, and 16 
particles. From table El we see an enormous growth of the standard shell-model diagonalization problem within the 
space. We calculated the number of M = states for He and B within the model space consisting of 4 major shells and 
estimated the number of basis states for C and O. Also indicated are similar estimates for seven major oscillator shells. 
The important lesson to learn from these numbers is that the model-space expansion becomes astronomical quite quickly. 



Table 9: Dimensions of the shell-model problem in four major oscillator shells and 7 major oscillator shells with M — 0. 



System 


4 shells 


7 shells 


"He 


4E4 


9E6 


8B 


4E8 


5E13 


12(j 


6E11 


4E19 


160 


3E14 


9E24 



Yet, because of the advent of radioactive nuclear beam accelerators, such as the proposed Rare Isotope Accelerator 
(RIA) in the U.S., we face the daunting task of moving beyond p-shell nuclei in ab initio calculations. We should therefore 
investigate several ways of approaching the nuclear many-body problem in order to successfully make the move into the 
RIA era. 
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One motivation for developing Auxiliary Field Monte Carlo (see section 0] above) for the shell model was to overcome 
the scaling problem. Other many-body methods that resum major classes of many-body diagrams have simply been 
neglected in nuclear science to this point. In this Section we will discuss the coupled-cluster technique which can be used 
to pursue nuclear many-body calculations to heavier systems beyond the p-shell. 

Coupled cluster theory originated in nuclear physics ^11201 around 1960. Early studies in the seventies probed 
ground-state properties in limited spaces with free nucleon-nucleon interactions available at the time. The subject was 
revisited only recently by Bishop et al. for further theoretical development, and by Mihaila and Heisenberg |36j . 

for coupled cluster calculations using realistic two- and three-nucleon bare interactions and expansions in the inverse 
particle-hole energy spacings. However, much of the impressive development in coupled cluster theory made in quantum 
chemistry in the last 15-20 years |22l 11261 11271 1128L 1129) still awaits applications to the nuclear many-body problem. 

Many solid theoretical reasons exist that motivate a pursuit of coupled-cluster methods. First of all, the method 
is fully microscopic and is capable of systematic and hierarchical improvements. Indeed, when one expands the cluster 
operator in coupled-cluster theory to all A particles in the system, one exactly produces the fully-correlated many-body 
wave function of the system. The only input that the method requires is the nucleon-nucleon interaction. The method 
may also be extended to higher-order interactions such as the three-nucleon interaction. Second, the method is size 
extensive which means that only linked diagrams appear in the computation of the energy (the expectation value of 
the Hamiltonian) and amplitude equations. As discussed in Ref. all shell model calculations that use particle-hole 
truncation schemes actually suffer from the inclusion of unconnected diagrams in computations of the energy. Third, 
coupled-cluster theory is also size consistent which means that the energy of two non-interacting fragments computed 
separately is the same as that computed for both fragments simultaneously. In chemistry, where the study of reactions is 
quite important, this is a crucial property not available in the interacting shell model (named configuration interaction 
in chemistry). Fourth, while the theory is not variational, the energy behaves as a variational quantity in most instances. 
Finally, from a computational point of view, the practical implementation of coupled cluster theory is amenable to parallel 
computing. 

We are in the process of applying quantum chemistry inspired coupled cluster methods [221 H'^OI 11261 11271 11281 11291 
I131L 11321 1133| to finite nuclei |7Tlll34| . We show one result from our current studies, namely the convergence of ^^O as a 
function of the model space in which we perform the calculations. 

The basic idea of coupled-cluster theory is that the correlated many-body wave function | 5*) may be obtained by 
application of a correlation operator, T, such that 

I =exp(-T) I $) , (179) 

where $ is a reference Slater determinant chosen as a convenient starting point. For example, we use the filled Gs state 
as the reference determinant for '*He. 
The correlation operator T is given by 

T = Ti+T2 + ---Ta, (180) 
and represent various n-particle-n-hole (np-nh) excitation amplitudes such as 

a{ef,i)ef 

T2 = \ J2 tt'alaUja,, (182) 

and higher-order terms for to Ta- We are currently exploring the coupled-cluster method at the Ti and T2 level. This 
is commonly referred to in the literature as Coupled-Cluster Singles and Doubles (CCSD). 
We compute the expectation of the ground-state energy from 

E^,,, = (*o I exp (-T) H exp (T) | ^-q) ■ (183) 

The Baker-Hausdorf relation may be used to rewrite the similarity transformation as 

exp (-T) H exp [T) = H + [H, Ti] + [H, T2] + ^ [[H, Ti] , Ti] + i [[H, T2] , T2] + [[H, Ti] , Ta] + • • • ■ (184) 

The expansion terminates exactly at four nested commutators when the Hamiltonian contains, at most, two-body terms, 
and at six-nested commutators when three-body potentials are present. We stress that this termination is exact, thus 
allowing for a derivation of exact expressions for the Ti {lp-\h) and T2 {2p-2h) amplitudes. The equations for amplitudes 
are found by left projection of excited Slater determinants so that 

= ($nexp(-r)i?exp(r)|$), 

= ($^'' I exp (-T) iJexp (T) | $) . (185) 
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The commutators also generate nonlinear terms within these expressions. To derive these equations is straightforward, 
but tedious, work [221 • While the resulting equations for the single and double excitation amplitudes appear quite lengthy, 
they are solvable through iterative techniques. 

The Hamiltonian may be written in a slightly more convenient form by explicitly calculating the expectation of the 
Hamiltonian in the reference state \ Eq = \ H \ This reference state is a single Slater determinant and 
represents, in this work, a doubly closed shell system. In this case, the Hamiltonian becomes 

H ^^fpq {alaq} + ^{pq\G\ rs) {fltaja^a,.} + Eq , (186) 

pq 

where the {} indicates normal ordering relative to the Fermi vacuum. The Fock operator is given by 

fpg = {p\K\q)+J2{P^\G\qt) . (187) 

i 

Using this form, and solving for the amphtudes, the energy of the system may then be calculated. The CCSD energy is 
(H) = EccsD = + \Y.('^\G\ ah)t'i^ + i ^(jj I G I ah)ttt] + E^ . (188) 

ia aibj aibj 

This equation is not restricted to the CCSD approximation. Since higher-order excitation operators such as T3 and T4 
cannot produce fully contracted terms with the two-body Hamiltonian, their contribution to the energy equation is zero. 
Higher-order excitation clusters can contribute indirectly to the energy through the equations used to determine the 
amplitudes. 

Because of the nonlinearity of the equations, one must have a good first guess for the np-nh amplitudes. For closed-shell 
nuclei, we use a M0ller-Plesset-Iike approach to generate the first guess for the iteration: 

- r) . ' 

t^}(X) - li^^i^^, (189) 

^ijab 

where Dia = fa — faa, and Dijab — fa + fjj — faa ~ fbb are single-particle energy denominators. Using this as a first guess, 
we can iterate the projection equations H185() until we find the converged amplitudes. Finally, we list the expressions for 
the Ti and T2 amplitudes. The Ti amplitude equations are given by 



c k kc kc kcd 

- 1 ^(fcz II ct)ti'{ -Y^hcm - Y.{ki II ci)titt + Y,{ka II cd)titf 

klc kc klc kcd 

- Y.^ki\\ cd)titttt + Y.^ki\\ cd)titt cd)tiiti -\Y.^ki\\ cd)t'fstf . (190) 



2 / J \ II " / Kl i <ry 

klcd klcd klcd klcd 

This equation is non-linear in the Ti amplitudes, and linear in the T2 amplitudes. 
The T2 amplitude equations are given by 

= (a6 1 1 zj) + ^ - fadt^) - ^ {h.t^t - Mt) 

c k 

+ \ Y.^kl II ^J)ttl + \ J2{ab II cd)t^, + Pi^mab)J2{kb \\ cj)C 

kl cd kc 

+ P(*j)^(a5 II cj)^^ - Piab)J2{kb \\ i])tl 

c k 

+ \p{^J)P{ab)Y,{kl II cd)t1^t'i^ + -JZ^kl\\ cd)t'^tl 

klcd klcd 

- ^^H) Y.^kl II cd)t^^tli - \P{^]) Y,{kl II cd)C't; 

klcd klcd 

+ lP{ab) Y.{kl\\ ij)tlt\ + \p{ij) Y^{ab\\ cd)i^i^ - P{ij)P{ab)Y^{kb II ic)tlt 



2fc 

''kl 



a^c 

2* \ — //_^\--" II "JZ-K-i ' 2' \"-''Z^\ — II —/"j-j ' \'J/- \"~//_^\-"- II --/"k'-j 
kl cd kc 
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+ P{ab) J2 fkctltf} + P{^]) E 

kc kc 

- P{ij) Y,{kl II ci)tlt^^ + P{ab) J2{ka \\ cd)tltf^ + P{ij)P{ab) Y,{^k \\ dc)tft% 

klc kcd kcd 

+ P{^J)P{ab)Y,{kl II icKt% + ip(zj)E(fc/ || cfc)t?i^f - ]^P{ab)Y^{kb \\ cd)tlt^f 

klc klc kcd 

- ip(zj)P(a6) J2{kb II cd)tXt'^ + lp{zj)P{ab) Y,{kl \\ cj)i!lif^ 

kcd klc 

- P(zj) Y,{kl\\ cd)tltX ~ P{ab) J2{kl\\ cd)tltttf^ + ^P{^J) Y,{kl\\ cd)tltpl\ 

klcd klcd klcd 

+ \p{ab)Y,{kl II cdnt^f^ + P{i])P{ab)Y,{kl II cd)t1t\tt^ 

klcd klcd 

+ \p{i])P{ab)Y,{kl\\cd)t^tltjt\ . (191) 



klcd 

The permutation operator P yields 



P{ij)f{ij) = f{ij) - f{ji) (192) 



The equations of the t2 amphtudes are nonhnear in both ti and T2 terms. While these equations appear quite lengthy, 
they are solvable through iterative techniques that we will discuss below. We note that the amplitude equations include 
terms that allow for 4p-4h excitations. Indeed, while we speak of doubles in terms of amplitudes, the class of diagrams 
involved in the th ----- - - . . - shell 

model with up to I s are linked 

and connected. Fc ^] 




CCSD Iteration 

Figure 26: The convergence of the ground-state energy as a function of the CCSD iterations in the ^^O system for various 
shell-model spaces. 



Shown in Fig.[5Blis the difference between the energy E{I) at a given iteration / and the final energy Ep as a function 
of iteration number for a calculation of ^^O in different oscillator spaces. Here we achieve convergence at the 10^^ level 
by 50 iterations in a model-space that includes seven major oscillator shells. 

Eq. I|189() represent the terms that one uses to compute the energy in second-order perturbation theory of the M0ller- 
Plesset type 23 . Shown in TableElare the energies obtained from the 0*'* E{Q) iteration and the final iteration E{F) as 
functions of increasing oscillator levels in the ^^O system. The difference between the converged CCSD energies and the 
initial 0*'* order energies increases as the basis space increases. The converged summation of the CCSD equations yields 
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Table 10: Comparisons of the 0*'' order energy E{0) and the converged CCSD results E{F) for ^^O as a function of 
increasing model model space. The results are also compared with many-body perturbation theory to second and third 
order, ^mbpt ^^'^ ^mbpT' respectively. All energies are in MeV. 



N (hcj) 


£;(o) 


Tp2nd 
''^MBPT 


pard 
-'^MBPT 


E{F) 


4 14 


-135.12 


-132.06 


-129.92 


-140.47 


5 14 


-124.79 


-124.84 


-121.52 


-127.79 


6 14 


-121.36 


-121.48 


-118.23 


-119.73 



approximately 10 MeV (or 0.6 MeV per particle) in extra binding at the Ucu minima. These findings are corroborated by 
those from many-body perturbation theory. It is therefore worth comparing these results with those from second-order 
and third-order many-body perturbation theory as well. These are labeled £^mbpt ^^"^ -^mbpt same table. The 

reader should notice that the zeroth iterations of the coupled-cluster schemes already includes corrections to the one- 
body amplitudes ti. However, the energy denominators used in the computation of the second-order diagrams of Fig.|5] 
(diagrams 2 and 3) have hole states determined by 

2 

3<F 

where F stands for the Fermi energy. We do not perform a self-consistent Brueckner-Hartree-Fock calculation however, as 
done by e.g., Gad and Miither .143. ■ The agreement with the zeroth order iteration and second-order perturbation theory 
is very good, especially for five and six major shells, as can be seen from Table However, for third-order perturbation 
theory one clearly sees fairly large differences compared with the coupled-cluster results. Typically, the relation between 
first- and second-order in perturbation theory for ^^O is given by a factor of '--^ 5 — 6. For e.g., = 5 and hoj — 14 MeV, 
we have —329.12 MeV from first order and —47.72 from second order. To third order we obtain a repulsive contribution 
of 3.32 MeV, to be contrasted with the almost 3 MeV of attraction given by higher-order terms in the coupled-cluster 
expansion. This indicates that many-body perturbation theory to third order is most likely not a converged result. 
An interesting feature to be noted from many-body perturbation theory calculations is that higher terms loose their 
importance as the size of the system is increased. For "^He the relation between first-order and second-order perturbation 
theory is given by a factor of ~ 3 — 4, depending on the value of hui. Calculations for ^''Ca not reported here indicate 
a relation of ~ 7 — 9 between first-order and second-order perturbation theory. This is somewhat expected since the 
G-matrix is smaller for larger systems, although the energy denominators become smaller. 

In the initial coupled-cluster study, we performed calculations of the ^^O ground state for up to seven major oscillator 
shells as a function of Tilo. Fig. |23 indicates the level of convergence of the energy per particle for A^ = 4, 5, 6, 7 shells. 
The experimental value resides at 7.98 MeV per particle. This calculation is practically converged. By seven oscillator 
shells, the fioj dependence becomes rather minimal and we find a ground-state binding energy of 7.52 MeV per particle in 
oxygen using the Idaho-A potential. Since the Coulomb interaction should give approximately 0.7 MeV/A of repulsion, 
and is not included in this calculation, we actually obtain approximately 6.90 MeV of nuclear binding in the 7 major shell 
calculation which is somewhat above the experimental value. We note that the entire procedure (G-matrix plus CCSD) 
tends to approach from below converged solutions. 

We also considered chemistry inspired noniterative triples corrections to the ground state energy. We performed this 
study in the model space consisting of four major oscillator shells. Table ITTI shows the total ground-state energy values 
obtained with the CCSD and one of the triples-correction approaches (labeled CR-CCSD(T) [12111123113311201 in the 
table). Slightly differing triples-corrections yield similar corrections to the CCSD energy. The coupled cluster methods 
recover the bulk of the correlation effects, producing the results of the SM-SDTQ, or better, quality. SM-SDTQ stands 
for the expensive shell-model (SM) diagonalization in a huge space spanned by the reference and all singly (S), doubly 
(D), triply (T), and quadruply (Q) excited determinants. To understand this result, we note that the CCSD Ti and 
T2 amplitudes are similar in order of magnitude. (For an oscillator basis, both Ti and T2 contribute to the first-order 
MBPT wave function.) Thus, the T1T2 disconnected triples are large, much larger than the T3 connected triples, and the 
difference between the SM-SDT (SM singles, doubles, and triples) and SM-SD energies is mostly due to TiT2.The small 
T3 effects, as estimated by CR-CCSD(T), are consistent with the SM diagonalization calculations. If the T3 corrections 
were large, we would observe a significant lowering of the CCSD energy, far below the SM-SDTQ result. Moreover, the 
CCSD and CR-CCSD(T) methods bring the nonnegligible higher-than-quadruple excitations, such as T1T2, JiTl, and 
T2 , which are not present in SM-SDTQ. It is, therefore, quite likely that the CR-CCSD(T) results are very close to the 
results of the exact diagonalization, which cannot be performed. 

These results indicate that the bulk of the correlation energy within a nucleus can be obtained by solving the CCSD 
equations. This gives us confidence that we should pursue this method in opened shell systems and to excited states. 
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Figure 27: Dependence of the ground-state energy of ^^O on huj as a function of increasing model space. 



Table 11: The ground-state energy of ^^O calculated using various coupled cluster methods and oscillator basis states. 



Method 


Energy 


CCSD 


-139.310 


CR-CCSD(T) 


-139.467 


SM-SD 


-131.887 


SM-SDT 


-135.489 


SM-SDTQ 


-138.387 



We have recently |134| performed excited state calculations on ''He using the EOMCCSD (equation of motion CCSD) 
method. For the excited states {"^k) and energies 

^(CCSD) > 0), we apply the EOMCCSD ("equation of motion 
CCSD") approximation (13211133] (equivalent to the response CCSD method fT37l), in which 

1*^) = i?(ccsD) exp(j.(ccsD))|^^_ (194) 

Here R'^^"""' ^ Rq + Ri + R2 is a sum of the reference (Ro), one-body (i?i), and two-body {R2) components obtained 
by diagonalizing /f(c;C!SD) -j^ ^j^g same space of singly and doubly excited determinants |$°) and l^^j") as used in the 
ground-state CCSD calculations. These calculations may also be corrected in a non-iterative fashion using the completely 
renormalized theory for excited states |128l 11291 11381 11361 1135| . The low-lying J = 1 state most likely results from the 
center-of-mass contamination which we have removed only from the ground state. The J = and J = 2 states calculated 
using EOMCCSD and CR-CCSD(T) are in excellent agreement with the exact results. 



Table 12: The excitation energies of '*He calculated using the oscillator basis states (in MeV). 



State 


EOMCCSD 


CR-CCSD(T) 


CISD 


Exact 


J=l 


11.791 


12.044 


17.515 


11.465 


J=0 


21.203 


21.489 


24.969 


21.569 


J=2 


22.435 


22.650 


24.966 


22.697 



Our experience thus far with the quantum chemistry inspired coupled cluster approximations to calculate the ground 
and excited states of the ^He and ^^O nuclei indicates that this will be a promising method for nuclear physics. By 
comparing coupled cluster results with the exact results obtained by diagonalizing the Hamiltonian in the same model 
space, we demonstrated that relatively inexpensive coupled cluster approximations recover the bulk of the nucleon corre- 
lation effects in ground- and excited-state nuclei. These results are a strong motivation to further develop coupled cluster 
methods for the nuclear many-body problem, so that accurate ab initio calculations for small- and medium-size nuclei 
become as routine as in molecular electronic structure calculations. 
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7 Three-body forces in shell-model studies 



An important feature of large scale shell-model calculations is that they allow one to probe the underlying many-body 
physics in a hitherto unprecedented way. As we have seen, the crucial starting point in all such shell-model calculations 
is the derivation of an effective interaction, be it either an approach based on a microscopic theory starting from the free 
NN interaction or a more phenomenologically determined interaction. In shell-model studies of e.g., the Sn isotopes, 
one may have up to 31 valence particles or holes interacting via e.g., an effective two-body interaction. The results of 
such calculations can therefore yield, when compared with the availiable body of experimental data, critical inputs to 
the underlying theory of the effective interaction. Thus, by going to the tin isotopes, in which the major neutron shell 
between neutron numbers 50 and 82 is being filled beyond the ^''^Sn closed shell core, we have the opportunity of testing 
the potential of large-scale shell-model calculations as well as the realiability of realistic effective interactions in systems 
with many valence particles. It should be noted that in many current shell-model calculations the effective interaction 
is frequently either parametrized or adjusted in order to optimize the fit to the data. As a matter of principle we shall 
refrain from making any such adjustments and stick to the interaction obtained by a rigorous calculation consistent with 
the many-body scheme chosen. Only then may one be able to assess the quality and reliability of the interaction obtained 
and the possible needs for improvement. It is our firm belief that one of the important aims behind many-body based 
derivations of effective interactions for the shell model is namely to provide a link between e.g., the free nucleon-nucleon 
interaction and properties of finite nuclei. Recalling our discussion in section ^ there are indications that three-body 
interactions, both real and effective ones, may be of importance. The Green's function Monte Carlo calculations of the 
Argonne-Urbana group 1111 IT^ and the recent no-core shell-model calculations O ll3Hj clearly indicate the need for 
three-body forces. Similar arguments, based on an analysis of Op, IsOd and IpOf nuclei by Zuker |14()| lend support to 
this picture for heavier nuclei as well. Real three-body forces are necessary in order to reproduce the saturation properties 
of nuclear matter as well, see for example Ref. |46| . 

Thus, with many valence nucleons present, such large-scale shell-model calculations may tell us how well e.g., an 
effective interaction which only includes two-body terms does in reproducing properties such as excitation spectra and 
binding energies. 

In general, excitation spectra for the chain of both odd and even isotopes from -'^^^Sn to ^^°Sn exhibit an excellent 
agreement with the data, see for example our discussion in section^ All these results are based on shell-model calculations 
starting with a two-body effective interaction determined from the recent charge dependent NN potentials of Machleidt 
and co-workers 64' '65' . These A^A^-interactions were renormalized for the given nuclear medium, using for example 
^^^Sn as closed shell core, through the introduction of the so-called reaction matrix G, which corresponds to solving the 
Lippmann-Schwinger equation for a finite nucleus. The G-matrix formed then the basis for a perturbative calculation 
of more complicated Feynman-Goldstone diagrams. However, if one studies Table ICT one sees that the binding energy 



Table 13: Binding energies for Sn isotopes. 





i^"Sn 


i^«Sn 


i^«Sn 


l^'^Sn 


^^^Sn 


^^^Sn 


"«Sn 


Exp 


-2.09 


-1.70 


-3.64 


-4.79 


-5.47 


-5.64 


-2.61 


1/(2) 


-2.24 


-2.63 


-4.60 


-6.99 


-9.39 


-11.77 


-18.58 


Mod. Shell Model 


-2.09 


-2.78 


-3.72 


-4.81 


-5.32 


-5.22 


-1.12 


+ ^(3) 




+0.009 


+0.029 


+0.096 


+0.213 


+0.394 


+0.998 


Trp. contr. 




+0.009 


+0.034 


+0.170 


+0.504 


+ 1.08 


+4.760 



relative to ^'^^Sn is clearly at askance with the data. The binding energy is defined as 

BE[^^'^-''Sn)] = BEi^'^^-'^Sn] - BE[^^^Sn] - n {BE[^^^Sn] - BE{[^^^Sn]) . (195) 

Experiment indicates a minimum around ^^"^Sn-^^^Sn and consequently a shell closure around ^^^Sn whereas theoretical 
binding energies increase linearly all the way down to ^^^Sn. 

The above is just an example of one of the problems which beset the theory of effective interactions for the shell 
model. In this case, as also pointed out by Zuker and co-workers |140| . one is not able to obtain simultaneously a 
good reproduction of both the excitation spectra and the binding energy. Similar problems have also been discussed in 
connection with large scale shell-model calculations of If Op shell nuclei. To give an example, effective interactions derived 
from two-body NN interactions which fit nucleon-nucleon scattering data, are not able to reproduce the well-known shell 
clousure in "'^Ca or the excitation spectra of ^^Ca and "^^Ca, see the discussion in Refs. |51 [TCl[Tin] . 

In this work we wish to address the discrepancy between theory and experiment shown in Table 1131 by including 
effective three-body forces in our shell-model calculations. The reason for this follows from the observation that the 
introduction of a global monopole correction of the form Wn{n — l)/2, with n being the number of valence particles and 
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W a quantity which is related to the energy centroids, can, when added to our theoretical binding energy, partly cure the 
discrepancy seen in Table IT^ Three-body forces typically yield repulsive corrections which scale as n(n — 1). In Table 
IIHI we adjusted the term W by simply requiring it to equal the difference in binding energy between the calculated and 
experimental values for ^'^°Sn. This resulted in W = 0.15 MeV. Adding such a global correction to even nuclei shown 
in Table [TSl results in the column labelled Mod. shell model. Clearly this improves the theoretical results in the correct 
direction and such a global modification of the matrix elements has no effect on the excitation spectra. 

These observations form thus the starting point for our investigation of three-body effective interactions. More 
explicitely, we aim at seeing whether three-body interactions may yield a more microscopic understanding of the above 
problem. The results presented employ effective three-body diagrams to second order in the G-matrix (see Polls et al. 
11421 for details) and have been shown at various conferences |141) . Examples of three-body diagrams are shown in 
Fig. I2H1 Diagram (a) is a third-order contribution not included in the Q-box. Diagram (b) is included. The inclusion of 




i 



(a) (b) 

Figure 28: Examples of unfolded three-body diagrams. 



such effective three-body diagrams do add a needed repulsion as seen from the column labelled V'^^'^ in Table [T51 However 
there is still large discrepancy between theory and experiment for the binding energy and even multiplying the number 
of triplets 

n\ 

number of triplets = — -- 

3!(n — 3)! 

with the three-body contribution from ^^^Sn yields too little attraction. Most likely the lack of repulsion can be ascribed 
to real three-body forces. Effective three-body forces for identical particles tend to be small due to the absence of the 
strong tensor channels in the intermediate states. This is also in line with recent works on the energy of pure neutron 
drops, where three-body clusters are included The reader should however note that our inclusion of three-body 
effective diagrams is incomplete. The summation of folded diagrams is based on only a two-body and three-body Q-box. 
We have not included a one-body Q-box. The latter introduces unlinked diagrams which need to be subtracted correctly 
when summing folded diagrams. To investigate whether the discrepancies seen for the binding energies are due to real 
three-body forces or not, we plan to include these in future studies, together with a three-body G-matrix and properly 
determined Q-box and folded diagrams. One needs thus to carefully distinguish between three-body forces and effective 
three-body interactions. 



8 Perspectives 

The study of exotic nuclei opens new challenges to nuclear physics. The challenges and the excitement arise because exotic 
nuclei will present new and radically different manifestations of nucleonic matter that occur near the bounds of nuclear 
existence, where the special features of weakly bound, quantal systems come into prominence. Furthermore, many of 
these nuclei are key to understanding matter production in the universe. Given that present and future nuclear structure 
research facilities will open significant territory into regions of medium-mass and heavier nuclei, it becomes important 
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to investigate theoretical methods that wiU aUow for a description of medium-mass nuclear systems. Such systems pose 
significant challenges to existing nuclear structure models, especially since many of these nuclei will be unstable and 
short-lived. How to deal with weakly bound systems and coupling to resonant states is an unsettled problem in nuclear 
spectroscopy. 

The aim of this work has been to expose both the weakness and strength of three many-body methods for studying 
nuclei heavier than ^^O. For lighter nuclei, both Green's function Monte Carlo [Sj IIOL IllL [T^ and no-core shell-model 
calculations |P 12 El 1 741 1139] offer true benchmarks for nuclear structure studies. The nuclear shell-model, combined 
with microscopically or phenomenologically derived effective interactions, has been extremely succesful in reproducing 
properties of nuclei employing Ohuj model spaces such as the IsOd the IpO/ shell [SlISS] and, as discussed in this work, 
medium heavy systems in the region of A = 100 to A = 132. However, if one wishes to study weakly bound systems 
or resonant states within this framework, the dimensionality of the problem soon exceeds present and most likely future 
hardware capabilities. 

For the heavier systems which are expected to be studied with facilities like RIA, many-body methods like the coupled 
cluster or parquet diagrams offer possibilities for extending microscopic ab-initio calculations to nuclei like ^"^Ca. Especially 
the coupled-cluster methods are very promising, since they allow one to study ground- and excited-state properties of 
nuclei with dimensionalities beyond the capability of present shell-model approaches. As demonstrated in Ref. |134| we 
show for the first time how to calculate excited states for a nucleus using coupled cluster methods from quantum chemistry. 
For the weakly bound nuclei to be produced by future low-energy nuclear structure facilities it is almost imperative to 
increase the degrees of freedom under study in order to reproduce basic properties of these systems. We are presently 
working on deriving effective interactions for weakly bound systems to be used in coupled cluster studies of these weakly 
bound nuclei. An extension of this project is to extract effective two-body interactions for open shell systems. These will 
in turn be used in shell calculations. 

We have based most of our analysis on effective interactions using two-body nucleon-nucleon interactions only. We feel 
this is important since techniques like the coupled cluster methods or parquet diagrams allow one to include a much larger 
class of many-body terms than done earlier. Eventual discrepancies with experiment such as the missing reproduction of 
e.g., the first excited 2+ state in a IpOf calculation of ''^Ca, can then be ascribed to eventual missing three-body forces, 
as indicated by the studies in Refs. |1(J[ 1111 IT^ I46L ITU I139| for light nuclei. The inclusion of real three-body interactions 
belongs to our future plans both for effective interactions for nuclei with A ~ 100 and larger as well as within the coupled 
cluster method and parquet diagrams for nuclei like ^°Ca. 
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